AVALIAÇÃO DO 
RISCO DE FALHA 


BASEADA NA MECÂNICA DA FRATURA 
E EFEITO DE ESCALA CONSIDERANDO 
A INCERTEZA DA CARACTERIZAÇÃO DE 


MATERIAIS CIMENTÍCIOS 


ORIENTADOR 
Francisco Evangelista Júnior / 


Y | 
APOIO 4 Sta 
E | / 
| ó 
/ 
14 


UnB CAPES 


AVALIAÇÃO DO 
RISCO DE FALHA 


BASEADA NA MECÂNICA DA FRATURA 
E EFEITO DE ESCALA CONSIDERANDO 
A INCERTEZA DA CARACTERIZAÇÃO DE 


MATERIAIS CIMENTÍCIOS 


Inés Damaris Munhoz Pena 


ORIENTADOR 
Francisco Evangelista Júnior |. 


Á 


y A 
Sa 

| | 

I 


APOIO 


Ina” CO> 


UnB CAPES Ano 2024 


A 


CIENCIAS EXATAS E DA TERRA E ENGENHARIAS 


Editora chefe 

Profº? Dr? Antonella Carvalho de Oliveira 
Editora executiva 
Natalia Oliveira 
Assistente editorial 
Flávia Roberta Barão 
Bibliotecária 

Janaina Ramos 
Projeto gráfico 

Ellen Andressa Kubisty 
Luiza Alves Batista 
Nataly Evilin Gayde 
Thamires Camili Gayde 
Imagens da capa 
iStock 

Edição de arte 

Luiza Alves Batista 


2024 by Atena Editora 

Copyright O Atena Editora 

Copyright do texto O 2024 Os autores 
Copyright da edição O 2024 Atena 
Editora 

Direitos para esta edição cedidos à 
Atena Editora pelos autores. 

Open access publication by Atena 
Editora 


Todo o conteúdo deste livro está licenciado sob uma Licença 
de Atribuição Creative Commons. Atribuição-Não-Comercial- 
NãoDerivativos 4.0 Internacional (CC BY-NC-ND 4.0). 


Goce 


O conteúdo do texto e seus dados em sua forma, correção e confiabilidade são de 
responsabilidade exclusiva dos autores, inclusive não representam necessariamente 
a posição oficial da Atena Editora. Permitido o download da obra e o compartilhamento 
desde que sejam atribuídos créditos aos autores, mas sem a possibilidade de alterá- 
la de nenhuma forma ou utilizá-la para fins comerciais. 


Todos os manuscritos foram previamente submetidos à avaliação cega pelos pares, 
membros do Conselho Editorial desta Editora, tendo sido aprovados para a publicação 
com base em critérios de neutralidade e imparcialidade acadêmica. 


A Atena Editora é comprometida em garantir a integridade editorial em todas as etapas 
do processo de publicação, evitando plágio, dados ou resultados fraudulentos e 
impedindo que interesses financeiros comprometam os padrões éticos da publicação. 
Situações suspeitas de má conduta científica serão investigadas sob o mais alto 
padrão de rigor acadêmico e ético. 


Conselho Editorial 

Ciências Exatas e da Terra e Engenharias 

Prof. Dr. Adélio Alcino Sampaio Castro Machado - Universidade do Porto 

Profº? Dr? Alana Maria Cerqueira de Oliveira - Instituto Federal do Acre 

Prof? Dr? Ana Grasielle Dionísio Corrêa - Universidade Presbiteriana Mackenzie 
Profº Drº Ana Paula Florêncio Aires - Universidade de Trás-os-Montes e Alto Douro 
Prof. Dr. Carlos Eduardo Sanches de Andrade - Universidade Federal de Goiás 
Profº Dr? Carmen Lúcia Voigt - Universidade Norte do Paraná 


A 


CIENCIAS EXATAS E DA TERRA E ENGENHARIAS 


Prof. Dr. Cleiseano Emanuel da Silva Paniagua - Colégio Militar Dr. José Aluísio da 
Silva Luz / Colégio Santa Cruz de Araguaina/TO 

Profº Drº Cristina Aledi Felsemburgh - Universidade Federal do Oeste do Pará 

Prof. Dr. Diogo Peixoto Cordova - Universidade Federal do Pampa, Campus Caçapava 
do Sul 

Prof. Dr. Douglas Gonçalves da Silva - Universidade Estadual do Sudoeste da Bahia 
Prof. Dr. Eloi Rufato Junior - Universidade Tecnológica Federal do Paraná 

Profº Dr? Érica de Melo Azevedo - Instituto Federal do Rio de Janeiro 

Prof. Dr. Fabrício Menezes Ramos - Instituto Federal do Pará 

Prof. Dr. Fabrício Moraes de Almeida - Universidade Federal de Rondônia 

Prof? Dr? Glécilla Colombelli de Souza Nunes - Universidade Estadual de Maringá 
Prof. Dr. Hauster Maximiler Campos de Paula - Universidade Federal de Viçosa 
Profº Drº lara Margolis Ribeiro - Universidade Federal de Pernambuco 

Prof? Drº Jéssica Barbosa da Silva do Nascimento - Universidade Estadual de Santa 
Cruz 

Prof? Drº Jéssica Verger Nardeli - Universidade Estadual Paulista Júlio de Mesquita 
Filho 

Prof. Dr. Juliano Bitencourt Campos - Universidade do Extremo Sul Catarinense 
Prof. Dr. Juliano Carlo Rufino de Freitas - Universidade Federal de Campina Grande 
Prof. Dr. Leonardo França da Silva - Universidade Federal de Viçosa 

Profº? Dr? Luciana do Nascimento Mendes - Instituto Federal de Educação, Ciência e 
Tecnologia do Rio Grande do Norte 

Prof. Dr. Marcelo Marques - Universidade Estadual de Maringá 

Prof. Dr. Marco Aurélio Kistemann Junior - Universidade Federal de Juiz de Fora 
Prof. Dr. Marcos Vinicius Winckler Caldeira - Universidade Federal do Espírito Santo 
Profº Dr? Maria laponeide Fernandes Macêdo - Universidade do Estado do Rio de 
Janeiro 

Profº Drº Maria José de Holanda Leite - Universidade Federal de Alagoas 

Profº Dr? Mariana Natale Fiorelli Fabiche - Universidade Estadual de Maringá 

Prof. Dr. Miguel Adriano Inácio - Instituto Nacional de Pesquisas Espaciais 

Prof. Dr. Milson dos Santos Barbosa - Universidade Tiradentes 

Profº Drº Natiéli Piovesan - Instituto Federal do Rio Grande do Norte 

Prof? Dr? Neiva Maria de Almeida - Universidade Federal da Paraíba 

Prof. Dr. Nilzo Ivo Ladwig - Universidade do Extremo Sul Catarinense 

Profº Dr? Priscila Natasha Kinas - Universidade do Estado de Santa Catarina 

Profº Drº? Priscila Tessmer Scaglioni - Universidade Federal de Pelotas 

Prof. Dr. Rafael Pacheco dos Santos - Universidade do Estado de Santa Catarina 
Prof. Dr. Ramiro Picoli Nippes - Universidade Estadual de Maringá 

Profº Drº Regina Célia da Silva Barros Allil - Universidade Federal do Rio de Janeiro 
Prof. Dr. Sidney Gonçalo de Lima - Universidade Federal do Piauí 

Prof. Dr. Takeshy Tachizawa - Faculdade de Campo Limpo Paulista 


A 


CIENCIAS EXATAS E DA TERRA E ENGENHARIAS 


Avaliação do risco de falha baseada na mecânica da 
fratura e efeito de escala considerando a incerteza 
da caracterização de materiais cimentícios 


Diagramação: 
Correção: 
Indexação: 
Revisão: 
Autores: 


Ellen Andressa Kubisty 

Maiara Ferreira 

Amanda Kelly da Costa Veiga 

Os autores 

Inés Damaris Muhioz Pefia 

Francisco Evangelista Júnior (Orientador) 


Dados Internacionais de Catalogação na Publicação (CIP) 


S598 Simonetti, Marcelo José 
Avaliação do risco de falha baseada na mecânica da fratura 
e efeito de escala considerando a incerteza da 
caracterização de materiais cimentícios / Inés Damaris 
Muhioz Pefia, Francisco Evangelista Júnior (Orientador). 
- Ponta Grossa - PR: Atena, 2024. 


Formato: PDF 

Requisitos de sistema: Adobe Acrobat Reader 
Modo de acesso: World Wide Web 

Inclui bibliografia 

ISBN 978-65-258-2649-3 

DOI: https://doi.org/10.22533/at.ed.493242006 


1. Mecânica da fratura. |. Pefia, Inés Damaris MuhÃioz. Il. 
Evangelista Júnior, Francisco (Orientador). III. Título. 


CDD 620.1362 


Elaborado por Bibliotecária Janaina Ramos - CRB-8/9166 


Atena Editora 

Ponta Grossa - Paraná - Brasil 
Telefone: +55 (42) 3323-5493 
www .atenaeditora.com.br 
contatoBatenaeditora.com.br 


A 


CIENCIAS EXATAS E DA TERRA E ENGENHARIAS 


DECLARAÇÃO DOS AUTORES 


Os autores desta obra: 1. Atestam não possuir qualquer interesse comercial que 
constitua um conflito de interesses em relação ao artigo científico publicado; 2. 
Declaram que participaram ativamente da construção dos respectivos manuscritos, 
preferencialmente na: a) Concepção do estudo, e/ou aquisição de dados, e/ou análise 
e interpretação de dados; b) Elaboração do artigo ou revisão com vistas a tornar o 
material intelectualmente relevante; c) Aprovação final do manuscrito para 
submissão.; 3. Certificam que o texto publicado está completamente isento de dados 
e/ou resultados fraudulentos; 4. Confirmam a citação e a referência correta de todos 
os dados e de interpretações de dados de outras pesquisas; 5. Reconhecem terem 
informado todas as fontes de financiamento recebidas para a consecução da 
pesquisa; 6. Autorizam a edição da obra, que incluem os registros de ficha 
catalográfica, ISBN, DOI e demais indexadores, projeto visual e criação de capa, 
diagramação de miolo, assim como lançamento e divulgação da mesma conforme 


critérios da Atena Editora. 


A 


CIENCIAS EXATAS E DA TERRA E ENGENHARIAS 


DECLARAÇÃO DA EDITORA 


A Atena Editora declara, para os devidos fins de direito, que: 1. A presente publicação 
constitui apenas transferência temporária dos direitos autorais, direito sobre a 
publicação, inclusive não constitui responsabilidade solidária na criação dos 
manuscritos publicados, nos termos previstos na Lei sobre direitos autorais (Lei 
9610/98), no art. 184 do Código Penal e no art. 927 do Código Civil; 2. Autoriza e 
incentiva os autores a assinarem contratos com repositórios institucionais, com fins 
exclusivos de divulgação da obra, desde que com o devido reconhecimento de autoria 
e edição e sem qualquer finalidade comercial; 3. Todos os e-book são open access, 
desta forma não os comercializa em seu site, sites parceiros, plataformas de e- 
commerce, ou qualquer outro meio virtual ou físico, portanto, está isenta de repasses 
de direitos autorais aos autores; 4. Todos os membros do conselho editorial são 
doutores e vinculados a instituições de ensino superior públicas, conforme 
recomendação da CAPES para obtenção do Qualis livro; 5. Não cede, comercializa ou 
autoriza a utilização dos nomes e e-mails dos autores, bem como nenhum outro dado 


dos mesmos, para qualquer finalidade que não o escopo da divulgação desta obra. 


Pá 


DEDICATÓRIA 


Salmo 135: 1 e 138:3 


“Agradeço a Deus, porque é infinitamente bom, 
porque a sua misericórdia dura para sempre, 
quando eu procurei a ele, ele respondeu e 
renovo minha força” 


Dedico este trabalho a minha Família que são 
fonte de amor e fortaleza e aos meus amigos. 


RESUMO 


Comumente, os sistemas estruturais são projetados com modelos determinísticos, 
desprezando a natureza incerta dos diversos parâmetros determinantes do 
comportamento estrutural. Além disso, boa parte dos modelos de falha estrutural 
utilizados nos projetos de dimensionamento não considera o racional mecânico, 
oriundo da mecânica da fratura na predição do colapso estrutural de sistemas e 
membros estruturais. Este trabalho estuda o risco estrutural das predições com 
modelos da Mecânica da Fratura Não-Linear Elástica para materiais e estruturas 
quasi-frágeis, especialmente para os materiais cimentícios com e sem adição de 
reciclados. Esta pesquisa implementa um algoritmo para realizar a análise de 
funções de estado limite não lineares que descrevem o estado limite de ruptura 
de vigas por propagação instável de trincas. Parâmetros de fratura de materiais 
cimentícios são extraídos de ensaios de laboratório como especificados nos 
métodos clássicos: efeito de escala de Bazant (EEB) e modelo de dois parâmetros 
(MDP). Os resultados quantificam a predição da falha estrutural e probabilidade 
de fraturamento instável incorporando a influência da variabilidade estatística 
dos parâmetros dos materiais obtidos em ensaios de laboratório. Os principais 
parâmetros de fratura investigados são a energia de fratura crítica (G,) e o tamanho 
equivalente da zona de processo (C). O trabalho faze comparações para medir 
a confiabilidade entre matérias tais concreto convencional (CONC) e concreto 
reciclado (FRAP). 


PALAVRAS-CHAVE: mecânica da fratura, análise de confiabilidade, efeito de 
escala. 


ABSTRACT 


Commontly, structural systems are designed with deterministic models, disregarding 
the uncertain nature of the various relevant parameters of structural behavior. In 
addition, many of the structural failure models used in design do not consider 
the fracture mechanics to predict the structural failure or collapse of systems 
and structural members. This dissertation studies the structural risk predictions 
considering models of Non-Linear Elastic Fracture for quasi-brittle materials and 
structures, especially for cement and recycled materials. A algorithm is proposed 
and implemented to perform the analysis of nonlinear limit state functions 
that describe the failure state for three-point bending beams for unstable crack 
propagation. Cementitious materials fracture parameters are extracted from 
laboratory tests as specified in the classical methods Bazant's size effect model 
(EEB) and two-parameter fracture model (MDP). The results quantify the prediction 
of structural failure and unstable fracture probability incorporating the influence of 
statistical variability of the parameters of the materials obtained in laboratory tests. 
The main investigated fracture parameters are the critical fracture energy (G,) and 
the equivalent size of the process zone (C). This work also compares the reliability 
index for beams of conventional concrete (CONC) and recycled concrete (FRAP). 


KEYWORDS: fracture mechanics, reliability analysis, size effect. 
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INTRODUÇÃO 


Materiais cimentícios são usados com muita frequência para a construção de membros 
e sistemas estruturais de pequeno e grande porte em todas as áreas da engenharia. O 
concreto, por exemplo, é o material de engenharia mais utilizado no mundo. Dessa forma, 
é importante conhecer as vantagens e desvantagens de suas principais propriedades e 
características, tais como resistência, durabilidade e custo. Ainda hoje, boa parte dos modelos 
de falha estrutural utilizados nos projetos de dimensionamento não considera a mecânica 
da fratura na predição do colapso estrutural de sistemas e membros estruturais compostos 
por materiais cimentícios. Apesar de relativamente bem documentado na literatura (Bazant 
e Planas, 1997; Shah et al., 1994; Li, 2011), o fraturamento do concreto é ainda um tópico 
desafiador devido a suas características especiais, tais como comportamento quasi-frágil, 
não aplicabilidade da mecânica da fratura linear elástica (MFLE), dependência do tamanho 
da estrutura e grande variabilidade estatística nas propriedades (Evangelista Jr et al., 2013; 
Evangelista Jr e Pera, 2014). 

O problema de falha por fraturamento está intimamente relacionado com análises 
estatísticas, através da teoria do “elo mais fraco” onde os estudos probabilísticos de Weibull 
criaram o modelo de mesmo nome, que explica a diminuição da resistência com o aumento 
do volume da estrutura. Os trabalhos clássicos de Bazant (1984) e Bazant e Kazemi (1990) 
formalizaram a lei do efeito de escala, ou efeito de tamanho, para estruturas de concreto, 
aqui chamada de Efeito de Escala de Bazant (EEB). O modelo é baseado em conceitos da 
mecânica da fratura não linear elástica (MFNLE), em modo | de fratura, levando em conta 
duas propriedades de fratura: energia inicial crítica de fratura (G) e tamanho equivalente da 
zona de processo (0). Aqui o efeito de escala na falha das estruturas prevê uma diminuição 
da resistência nominal, a medida que o tamanho da estrutura aumenta. Paralelamente, 
Jenq e Shah (1985) utilizaram os mesmos conceitos da MFNLE para propor o Modelo de 
Dois Parâmetros (MDP), que também é capaz de prever o efeito de escala em estruturas de 
concreto por meio do G,e do deslocamento crítico da ponta da trinca (CTOD,). 

Os dois modelos EEB e MDP vêm sendo os principais modelos de previsão de 
crescimento de trincas em estruturas de concreto, inclusive com métodos RILEM 
diferenciados para o cálculo de seus parâmetros (RILEM 1990a; RILEM 1990b). Desde 
então, estes modelos e seus respectivos procedimentos de ensaio tem criado muita 
controvérsia, em termos de eficiência e facilidade da extração dos parâmetros de fratura 
em laboratório e acurácia de suas previsões (Mehta e Monteiro, 2012; Mindess et al., 20083; 
Bazant; Planas, 1997; Shah et al., 1994; Li, 2011). 

O estudo ou avaliação do risco de falha baseia-se de forma geral em análises de 
confiabilidade, para medir a probabilidade de ocorrência em que um sistema ou membro 
estrutural irá desempenhar sua função no limite da falha estrutural. A confiabilidade 
é expressa quantitativamente por um índice de confiabilidade B, que se relaciona 


Introdução 


inversamente com a probabilidade de falha (P) da estrutura. A análise de confiabilidade 
considera as possíveis incertezas apresentadas nas variáveis do problema, para predizer 
a probabilidade da falha do sistema. Por meio da teoria de probabilidade e de uma função 
que relaciona estas variáveis (função de estado limite) têm-se as bases para estimar a 
confiabilidade através de algoritmos mais simples, tais como o primeira ordem segundo 
momento (First-Order Second Moment, FOSM), até métodos mais elaborados como o 
método de confiabilidade de primeira ordem (First-Order Reliability Method, FORM) (Nowak 
e Collins, 2013; Melchers (1999). Dependendo da forma da função do estado limite e do 
tipo de distribuição de probabilidade das variáveis envolvidas, soluções analíticas podem 
ser encontradas. Para casos mais complexos, principalmente envolvendo não linearidade 
da função no limite de falha, métodos numéricos mais sofisticados devem ser empregados 
para o cálculo do índice B. Também poderão ser analisados os parâmetros Ce G,e 
estabelecer uma comparação por meio do Método de Monte Carlo (MMC). 


1.1 JUSTIFICATIVA 


Comumente, os sistemas estruturais são projetados considerando modelos 
determinísticos, que negligenciam a natureza incerta dos diversos parâmetros determinantes 
do comportamento estrutural. No caso da predição probabilística da falha da estrutura por 
propagação instável de trincas, é de vital importância as variáveis de carregamento, de 
geometria e propriedades dos matérias. Além disso, a escolha do modelo a ser utilizado 
bem como as propriedades de fratura associados a ele assumem fundamental importância 
na predição probabilística da falha. Entretanto, nenhum trabalho da literatura investiga os 
efeitos dos dois principais modelos para materiais quasi-frágeis (EEB e MDP), considerando 
métodos probabilísticos para avaliar as incertezas na determinação dos parâmetros de 
fratura e como elas impactam a predição dos índices de segurança (e análise de risco) de 
membros estruturais. Também é de fundamental importância o estudo de como o efeito de 
escala influencia na confiabilidade estrutural considerando a falha por propagação instável 
de fissuras que levam ao colapso. 


1.2 OBJETIVOS 


O objetivo geral do presente trabalho é implementar o algoritmo para análise de 
funções de estado limite não lineares que descrevam a ruptura de vigas por propagação 
instável de trincas. Com isso, objetiva-se a quantificação dos efeitos das diferentes 
metodologias, e as incertezas dos parâmetros, de dois métodos clássicos de fratura de 
materiais cimentícios (EEB e MDP) nas predições de confiabilidade estrutural. 

Os objetivos específicos desta dissertação são: 
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*— implementar algoritmos eficientes para análises de estados limites não lineares 
por meio de versões avançadas do FORM. 


* avaliar o risco de falha de uma estrutura por meio de um algoritmo de confiabi- 
lidade e quantificar o índice de segurança considerando conceitos de mecânica 
de fratura; 


* estudar a relação entre índices de confiabilidade e modelo de efeito de escala 
em vigas em flexão em três pontos de concreto; 


* verificar a influência da consideração de variáveis aleatórias de carregamento, 
geometria e propriedades dos materiais no cálculo do índice de confiabilidade; 


* considerar os estados limites da resistência de com vários tamanhos e geo- 
metrias incorporando os modelos de EEB e MDP e como se comparam com a 
mecânica da fratura linear elástica; 


* determinar os índices de confiabilidade para vigas de concreto convencional 
(CONC) e concreto reciclado (FRAP); 


* quantificar a influência da variabilidade estatística dos parâmetros dos materiais 
obtidos em ensaios de laboratório de concreto convencional e reciclado na pre- 
dição da falha estrutural; 


* | realizar simulação de MMC para determinar os tipos de distribuições que apre- 
sentam os parâmetros mais importantes relacionados com os modelos de EEB 
e MDP para o CONC. 


1.3 CONTRIBUIÇÕES 


A pesquisa contribui para o estado da arte da análise de confiabilidade da falha de 
estruturas por propagação instável de trincas utilizando modelos da mecânica da fratura 
não linear elástica para materiais quasi-frágeis tais como materiais cimentícios. O trabalho 
contribui para uma análise mais profunda sobre os dois métodos mais utilizados para 
fratura em concreto (EEB e MDP) e de como suas respectivas metodologias influenciam 
na predição probabilística de falha estrutural. Além disso, os resultados contribuem para a 
análise da dependência dos índices de confiabilidade com o efeito de escala e propriedades 
fundamentais do material determinados por modelos de fratura. 

A metodologia empregada permite fazer comparações dos índices de segurança B 
entre dois materiais de muita utilização na construção atual, tais como concreto comum 
CONC e concreto reciclado FRAP, determinando qual pode representar maiores níveis 
de segurança em vigas. Os resultados obtidos poderão auxiliar na elaboração de projetos 
baseados em mecânica da fratura, pois determinarão as melhores combinações de 
geometria e propriedades de materiais de uma viga. 
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1.4 RESUMO DA METODOLOGIA 


O presente trabalho inicialmente implementa um algoritmo para análise de 
funções de estado limite não lineares que descrevam o estado limite de ruptura de vigas 
por propagação instável de trincas. Os parâmetros de fratura de materiais cimentícios 
convencionais e reciclados foram determinados de ensaios de laboratório como 
especificados nos métodos RILEM para os dois modelos clássicos EEB e MDP. A partir 
destes dados, as análises de confiabilidade consideram funções do estado limite de ruptura 
e as incertezas de variáveis geométricas, de carregamento e de propriedades de fratura 
de dois materiais: concreto convencional (CONC) e concreto reciclado (FRAP) em vigas de 
flexão em três pontos. A confiabilidade da resistência nominal de vigas de vários tamanhos 
são analisadas. Finalmente, os tipos de distribuição de probabilidade da resistência das 
vigas são determinadas por meio de simulações do método de Monte Carlo para que 
apresentam os parâmetros mais importantes relacionados com os modelos de EEB e MDP 
para o concreto. 


1.5 ORGANIZAÇÃO DA DISSERTAÇÃO 
O trabalho está organizado nos seguintes capítulos: 


* CAPITULO 2: Revisão bibliográfica de mecânica da fratura, conceitos básicos 
de estatística, confiabilidade e métodos de confiabilidade. 


*— CAPITULO 3: Neste capitulo, é descrita a metodologia utilizada para a avaliação 
de confiabilidade realizada para a flexão de vigas em três pontos (FTP) e ela- 
boradas com concreto convencional (CONC) e com concreto reciclado (FRAP). 


* CAPITULO 4: Neste capitulo apresentam-se os resultados das análises de con- 
fiabilidade realizados para o concreto convencional CONC e para o concreto 
reciclado FRAP. 


* CAPITULO 5: Neste capítulo apresenta as funções de densidade de probabili- 
dade da resistência de vigas em FTP usando a técnica de simulação do método 
de Monte Carlo (MMC). 


«CAPITULO 6: Apresenta as conclusões, aportes e sugestões para trabalhos 
futuros relacionados com a pesquisa desenvolvida. 
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REVISÃO BIBLIOGRÁFICA 


Neste capitulo apresenta-se a revisão bibliográfica para Mecânica da fratura 
linear elástica e não linear elástica da seção 2.1. até a seção 2.3. e para os conceitos de 
confiabilidade da seção 2.4. até a 2.7. 


2.1 MECÂNICA DA FRATURA 


A mecânica da fratura descreve a cinética e cinemática de corpos contendo trincas 
e a consequente propagação destas. O processo de faturamento, ilustrado na Figura 2.1, 
pode ser classificado em três tipos ou modos; 


Modo | Modo Il Modo III 
Por abertura Por cisalhamento no plano Por cisalhamento fora do plano 


Figura 2.1. Modos de fratura 


A lei do efeito de escala para estruturas de concreto para o modelo EEB é usada 
para conceitos da mecânica da fratura não linear elástica (MFNLE), em modo | de fratura, 
levando em conta as propriedades de fratura (G) e (C). Os mesmos conceitos da MFNLE são 
usados para propor o MDP, que também prevê o efeito de escala em estruturas de concreto 
por meio do G, e do deslocamento (CTOD,). Por outra parte, o conceito de mecânica da 
fratura linear elástica (MFLE), que assume uma deformação plástica na ponta de fissura, é 
desenvolvida por Griffith (1920) e também segue a lei de EEB para o parâmetro de fator de 
intensidade de tensões (FIT). 

As seções a seguir descrevem as abordagens para a descrição do processo de 
iniciação e propagação de trincas utilizadas neste trabalho. 


2.2 MECÂNICA DA FRATURA LINEAR ELÁSTICA (MFLE) 


Os critérios usados pela mecânica da fratura linear elástica (MFLE) são baseados 
em uma propriedade independente, que descreve a resistência ao fraturamento (ou 
tenacidade) em materiais frágeis: o fator de intensidade de tensões (FIT ou K) ou taxa 
de liberação de energia (TLE ou G,) se usadas as formulação de tensões ou de energia, 


respectivamente. 
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2.2.1 Fator de intensidade de tensão (FIT ou K) 


As imperfeições geométricas na superfície de um corpo gera uma distribuição 
diferente de tensões. Irwin (1958) indica que o FIT depende da geometria e da tensão 
induzida (0) pela carga atuante no corpo sólido. Para este trabalho em particular se 
considerado somente o modo | de fraturamento (por abertura). Neste caso, o FIT (K) pode 


ser expresso como: 
a 
Kj=ovtaf (5) (2.1) 


Em que a, é o comprimento inicial da fissura, D é a dimensão do corpo na direção 
da fissura; e f(a,/D) uma função que relaciona a, e D denominada como fator geométrico, 
que depende da geometria de cada corpo em questão. Para uma placa submetida a uma 
tensão uniaxial, como ilustra a Figura 2.2 e a Equação (2.2): 


ON 


Figura 2.2. Placa submetida a uma tensão uniaxial 


f (5) = |sec (55) (2.2) 


Para tensão uniaxial com eixo simples de fissura ilustrada na Figura 2.3 e na 
Equação (2.3) tem-se: 
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ON 


Figura 2.3. tensão uniaxial com eixo simples de fissura 


2 3 4 


(5) = 112-0.2315+ 10.55 (5) 21570 (5) + 30.39 (5) (2.3) 
“- 199-(5)(1-3) ) fp as - 393 (5) +2.70(5 | foi 


3 
vr(1+25)(1 5) 
A Equação (2.4) corresponde a uma viga de flexão em três pontos ilustrada na 
Figura 2.4: 


Figura 2.4. viga de flexão em três pontos 


É importante mencionar que Tada et al. (2000) apresenta r(5) para diversas 
geometrias e carregamentos. 
A Equação (2.1) pode ser expressa em termos da tensão máxima induzida 
(resistência do corpo) como: 
Kic 


o 


0, = 


(2.5) 
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Em que K,, é denominado de fator de intensidade crítico ou tenacidade à fratura 
(no modo |) que é considerado uma propriedade fundamental do material. Note que este 
fator é o valor limite (crítico) para que a trinca de tamanho a se propague instavelmente. 
Consequentemente, a tensão induzida o, é a máxima tensão suportada pelo corpo 
(resistência) para que occorra a propagação instável. 


2.2.2 Taxa de liberação de Energia de fratura 


Griffith (1921) estabelece a descrição da MFLE por meio do conceito de energia de 
propagação de fissuras, G,, que é também conhecida como a taxa de liberação de energia. 
Esta é também função da geometria da estrutura e da tensão induzida pelo carregamento. 
Para este estudo denomina-se G, como para modo | de fratura. Além disso, G, é também 


relacionado com K, para materiais elásticos, por meio da expressão: 


K, 
G== (2.6) 
E 


Em que E = E para o caso de estado plano de tensão; e E = E/(1 —- v?) para o 
estado plano de deformação; E é o módulo de elasticidade; e v é o coeficiente de Poisson. 
Do mesmo modo, uma propriedade do material denominada energia de fratura (G,= Go); 
que é a energia G, crítica para a propagação instável, e é relacionada com o fator de 
intensidade de tensão crítico por: 


Kic = (G,E)!? (2.7) 


O K, é assumido instável para o caso de estado plano de tensão onde E = E. 
Substituindo a Equação (2.1) temos: 


1/2 
gem (2.8) 


Vas (5) 

em que o, é a tensão máxima induzida no corpo para que occorra a propagação 
instável da trinca. 

Apesar de largamente aplicada, a MFLE só é adequada para materiais elásticos e 
que não exibam nenhum, ou desprezível, comportamento não linear na ponta da fissura. A 
maioria dos materiais da engenharia, exibe comportamento não linear na ponta da fissura, 
denominada zona de processo, em condições normais de operação. Nestes materiais a 
MFLE não pode ser aplicada, e teorias mais sofisticadas devem ser utilizadas. No caso de 
materiais cimentícios, a zona de processo na frente da trinca não é desprezível, e faz com 
que esses materiais apresentem comportamento quasi-frágil. 
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2.3 MECÂNICA DA FRATURA NÃO LINEAR ELÁSTICA (MFNLE) 


O concepto de trinca equivalente é abordado por dois modelos representativos da 
mecânica da fratura que são o (MDP) e (EEB). Estes modelos de trincas equivalentes, em 
lugar do comprimento original da trinca utilizam um comprimento de trinca equivalente que 
pode prescrever-se explicitamente como o parâmetro (C), e Wells (1961) evidencio que as 
faces da trinca movimentavam-se com deformação plástica o que permitiu-lhe desenvolver 
um parâmetro que conhecido como o deslocamento da apertura do frente de trinca (CTOD). 

A equivalência é possível devido a que o comprimento da trinca é conhecida e 
estável antes que o concreto seja submetido ao carregamento máximo, e o comprimento 
critico equivalente para a falha do concreto é diferente ao comprimento inicial da trinca. 
Para predizer a falha do concreto é necessário conhecer a carga crítica e o comprimento 


de trinca correspondentes 


2.3.1 Modelo do Efeito de Escala de Bazant (EEB) 


O modelo do Efeito de Escala de Bazant (EEB) analisa o comportamento não linear 
da fratura em uma estrutura de concreto, com uma força ou tensão de carregamento 
máximo, em função do tamanho da estrutura. A derivação detalhada do modelo pode ser 
encontrada em (Bazant, 1984; Bazant e Kazemi, 1990; Bazant e Planas, 1997). Entretanto 
a principal equação do EEB relacionando a resistência nominal de um membro estrutural, 
O, com as propriedades de fratura do material, dimensões e geometria da estrutura taxa de 
liberação de energia de fratura G,e o comprimento crítica da região de fratura C é: 

1 
GE 2 
= (Dre) 28) 

Em que C é uma propriedade do material que representa o comprimento equivalente 

da ZPF do material; D é uma dimensão característica (tamanho) da estrutura,a = = a é 


um comprimento de trinca na estrutura, e g(a) é um fator geométrico dado por: 
g(a) = Malf(a)C,]? (2.10) 
Em que C, é uma constante que relaciona o, e D que depende das dimensões e tipo 
de geometria; e f(a.) é a função geométrica da MFLE. É importante salientar que a Equação 


(2.9), devido á forma de E, é somente válida para o estado plano de tensão. 
Onde g'(a) faze referencia á derivada da função g(o.) descrita na Equação (2.10). 
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Log (Resistência) náo 
Esqemlsalosds sds N Critério de 
resistência 


Amostras lab. 
testes 


! tamanho pequeno 


Tamanho grande 
Mais estruturas 


Deflexão relativa Log (Escala/Tamanho) 


(a) curvas do ensaio (b) curvas de resistência com efeito de escala 


Figura 2.5. Modelo do efeito de escala de Bazant (EEB) com : a) curvas do ensaio de resistência 
nominal vs. deflexão relativa, e (b) escala logarítmica da resistência nominal vs. escala / tamanho . 
Adaptado de ACI comitê 446 1992) 


O (EEB) é definida pela comparação de estruturas de geometria similar de diferentes 
tamanhos. É importante definir que o termo o, representa o valor da resistência nominal 
no máximo valor do carregamento antes da falha instável carregamento (último) P. A 
Figura 2.5 para a estruturas geometricamente similares de diferentes tamanhos descreve 
o comportamento típico das curvas de esforço vs. deflexão, como é ilustrado a falha 
ocorre perto do pico, com os incrementos de tamanho. Este efeito pode predizer a fratura 
mecânica a partir do facto que para estruturas mais grandes tem-se uma maior energia de 
deformação que incrementa a propagação da zona de falha. 

Em que o, para a lei de efeito escala pode ser expressa como: 


Bfc 


"VIAD/D, (2.11) 


GE (1 
g'(0)C 


Bf' = | (2.12) 


g'(a) 
gta) “ 


Dj= (2.13) 

Em que ff é um esforço de tensão do material usado apenas para fins dimensionais, 
B é uma constante adimensional, e D, é uma constante com a dimensão de comprimento. 
Be D, dependem de das propriedades das propriedades de fratura do material e da fator 
geometria da viga. 
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On artantteart GR 


AAA AAAAAS! AAA A AAA AMRS 
Y 


Es 
Wi nvy RR 
(a) trinca blunt (b) trinca sharp 


Figura 2.6. Idealização da zona de processo para a) trinca blunt b) trinca sharp. Adaptado de ACI 
comitê 446 1992. 


Na Figura 2.6 a zona de processo de fratura considerada na ponta de uma trinca 
contínua, de extensão finita pode, por aproximação, ser considerada constante. Para a 
trinca sharp da Figura 2.6b, tensão normal transversal diminui gradualmente a partir de f/ 
até O ao longo do comprimento da zona de processo. Por causa da presença desta zona da 
trinca elasticamente equivalente, gera-se a liberação de energia de deformação do material 
que produz um incremento no comprimento efetivo de trinca a,, e que é assumido por o 
comprimento C, considerado como uma constante do material. 

O EEB também é usado apara a determinação das propriedades de fratura do 
material, G,e C. A metodologia de ensaio é normatizada em uma norma (RILEM 1990a) 
para vigas com entalhe central, a,, submetidas à flexão em três pontos FTP como ilustrado 
pela Figura 2.7. As vigas têm altura D, vão Se espessura t Basicamente o ensaio consiste 
em medir os máximos carregamentos, P, de vigas de três alturas/tamanhos (D), mas 


geometricamente similares (mesmos ). 


E Em 


Figura 2.7. Configurações do ensaio do EEB para vigas em flexão em três pontos. Adaptado de Bazant 
et al. (1997) 
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O carregamento máximo final (PP) pode ser obtido tomando o peso do corpo de 
prova em conta. 
PP=B+5LW (2.14) 
j 
Em que P, é o carregamento máximo medido no ensaio, S, o vão, L, comprimento 
total e WW,o peso da viga j. 


2 
g= (5) (2.15) 


J 
Fazendo X, = D, para j = 1, 2, 3,..,n, podemos estabelecer uma regressão linear 
definida como Y = Ax + B,. Os valores da inclinação A e o ponto de intercepção B, são: 


A = Bis) 


2 (2.16) 
Tal; e x) 


Bo=Y-AX (2.17) 


. 1% 
X= 4) (2.18) 
j=1 
e 
Ro 
F="50) (2.19) 
j=1 


Para o caso das vigas da Figura 2.7, o valor de fw,), e consequentemente g(a,), 
depende da razão S/D. O fator geométrico f(«,) para vigas com S/D = 4 é: 


ay 199-(58) (1-0) [pas - 3.93 (58) +2.70(99) | 
f(go = 38) = ADA ARA] 
m(1+27)(1-D) 
Os valores de G, podem ser determinados por: 


— g(ao) 
fo EA 


(2.20) 


(2.21) 


O valor do comprimento equivalente da ZPF, C, é por sua vez obtido por: 


— 9(00)Bo 


= FICAVI (2.22) 


O carregamento máximo medido no ensaio (P) se relaciona com a resistência 


nominal da viga em flexão em três pontos (FTP), com a altura D pela relação: 
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Or; = o (2.23) 
Em que: 
35 
Cn = 30 (2.24) 


A Equação do EEB pode também ser escrita na forma 


di 
Ch ( GE E (2.25) 


o = E E 
á tD Dg(ao) + Cg'(ao) 


Onde P, tem o mesmo sentido que P, é o máximo carregamento suportado pela 
viga, ou seja, sua resistência nominal. Analogamente à MFLE também pode fornecer uma 
relação para a resistência nominal adaptando-se a Equação (2.8): 


1 
GB (GEN E 
a o! (E e8) 


2.3.2 Modelo de dois parâmetros (MDP) 


Jeng e Shah (1985) adaptaram o trabalho de Wells (1961) seu modelo considera 
que a falha estrutural ocorre, quando (i) o FIT atinge o valor crítico K, e (ii) a abertura da 
ponta da fissura, aqui chamada de CTOD (crack tip opening displacement), atinge um valor 
crítico CTOD,.. Ou seja: 


K = Ki (2.27) 


CTOD = CTOD,, 


É assumido que, em um corpo sólido com entalhe a, estes dois valores críticos 
são atingidos para um comprimento de trinca equivalente a, = a, + Aa, em que Aa é 
um comprimento virtual. Pelo fato de usar duas propriedades de fratura (K, e CTOD,) 
este modelo é chamado de modelo de dois parâmetros (MDP). O método usa a curva 
de carregamento (P) e a abertura da boca da fissura, aqui chamada de CMOD (crack 
mouth opening displacement). Uma norma (RILEM 1990b) foi desenvolvida para materiais 
cimentícios utilizando ensaios de flexão em três pontos como ilustrado na Figura 2.8.a. A 
determinação dos parâmetros K e CTOD,, é feita a partir da medida das inclinações das 
curvas iniciais de carregamento (C) e descarregamento (C ) utilizando as Equações da 
MFLE. 
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Ho 


CcMOD 


(a) Configuração da viga 


Descarga no pico de carga 


CMOD 


CMOD$ CMODE 


(b) Curvas de carregamento e descarregamento 


Figura 2.8. Ensaio para o MDP: (a) configuração da viga de flexão em três pontos; (b) curvas de 
carregamento e descarregamentos e respectivas inclinações. 


Na derivação das equações do MDP, o CMOD,, é dividido em duas partes: uma 


parcela elástica CMODS; e uma parcela plástica CMODP, de tal foma que: 
CMOD,, = CMODe + CMODP? (2.28) 


Os módulos de elasticidade do material no processo de carregamento e 


descarregamento, E, e E,, respectivamente, podem sere determinados a partir das 
inclinações C, C, (Figura 2.8. b) e as dimensões do corpo de prova nas relações da MFLE: 
BOB (2.29) 


= cD2ç fo(%o) 


E; 


Era (2.30) 


“ CD2t 


E fo (ac) 
Em que D é a altura da viga, t é a espessura da viga e: 


ao +H, 
o+ lo (2.31) 


ESA 


e 
Ac + Ho 
(2.32) 


(ac) = Ho 
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Em que H é a altura da haste de suporte. 
O fator geométrico f(a.) e f(a,) é o mesmo nas duas expressões e é definido como: 


0.66 


ray (2.33) 


fola) = 0.76 — 2.28 « + 3.784? — 2.040º + 


Igualando as Equações (2.29) e (2.30), tem-se o comprimento de fissura equivalente, 


a Culog2(%o) (2.34) 
Cigo(aç) 
Agora, os valores de K,.e CMODS são calculados por meio das equações da MFLE 


considerando o comprimento de trinca efetiva a e a tensão máxima o: 
Kic = OTA AC) (2.35) 


Em que f(«,) é uma função adimensional geométrica. 


40, 


a 
CMODE = E “ fila) (2.36) 


Em que f(«.,) é uma função adimensional geométrica e E é o módulo de elasticidade. 
As Equações (2.35) e (2.36) podem serem formuladas em termos do carregamento 
máximo, P,, como: 


Sytacfila 
Kic = 3(B. + a (2.37) 


Em que W, = “2º com W, sendo peso da viga. 


6(P. + 0.5Wp)Sacfo(ac) 


2.38 
ED2t 


CMODE = 


O valor do deslocamento crítico da ponta da trinca CTOD, é determinado por; 


CTOD, = CMOD* fa, y,) (2.39) 


Em que y, = y/a,. 
As funções geométricas, f, e f, para vigas submetidas à flexão em três pontos são 
fornecidas por Murakami et al. (1987) como: 
1.99 — (e )(1 — «.)[2,15 — 3.93(00) + 2.70(0,)2] 
VTC +20) - a 


fila) = (2.40) 
1.73 — 0.85(a,) + 31.2(0,)? — 46.3(0,)? + 2.70(0,)* 


a E (2.41) 
— a) 


fla) = 
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fslac vo) = (1 — yo)? + (1.081 — 1.14900) (yo — vo?) (2.42) 


De acordo com Shah et al. (1994), o comprimento equivalente da ZP () utilizado no 
EEB pode ser relacionado com o CTOD,, por meio da expressão da viga de flexão em três 
pontos FTP: 


0.061 CTOD2E >» CTODZEa, CTOD2E 
C="— "05440 + [029249 ++ 00038(—— 2 (2.43) 
f f f 


Em que G, é calculado por meio da MFLE com a Equação (2.7) e K,. determinado 
pela Equação (2.35) ou (2.37) 
Então, a Equação (2.25) da predição do EEB pode ser usada com os valores de G, 


e C calculados a partir do K, e do CTOD, determinados pelo MDP. 


2.4 ANÁLISE DE CONFIABILIDADE 


Na engenheira o projeto de dimensionamento por exemplo, deve proporcionar aos 
elementos e aos sistemas estruturais critérios de segurança e de serviço da estrutura baixo 
à demanda de durabilidade. Além disso boa parte dos parâmetros relacionados com o 
design sempre apresentam numerosas incertezas, pelo que é necessário quantificar estas 
incertezas por meio de métodos de confiabilidade, que medem estas incertezas através 
de índices de confiabilidade e determinam os níveis de segurança, identificando os 
porcentagens de probabilidade que tem o sistema ou elemento de falhar. 

Alguns de os autores da literatura probabilística que estudam diferentes métodos 
de confiabilidade são Haldar e Mahadevan (2000), Hasofer e Lind (1974), Melchers (1987), 
Choi, Canfield e Grandhi, (2005) e Nowak e Collins, 2018. 

Haldar e Mahadevan (2000) fazem o desenvolvimento das etapas principais para 
quantificação de incertezas, em casos que as fontes de incerteza são incógnitas, sendo 
necessária informação adicional. Por exemplo, para avaliar a incerteza de qual é a máxima 
velocidade do vento de uma cidade em um ano, precisa-se reunir os dados dos maiores 
valores encontrados em anos anteriores, que podem ser reunidos em n grupos de anos, 
para n amostras. Inicialmente, são necessárias informações estatísticas que podem ser 
obtidas do quadro mostrado na Figura 2.9. A informação reunida é chamada espaço de 
características de aleatoriedade e pode ser representado em um histograma de frequência. 
O diagrama de frequência pode ser obtido de um jeito mais simples, aplicando a teoria de 
Função de Densidade de Probabilidade (FDP). 
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Mundo Real 


Informações importantes 
Espaço amostral 


Representação matemática de 
incertezas 


Densidade ou função de distribuição Estatísticas ou parâmetros de estimativa 


Definir os critérios de 
desempenho 


y 
Avaliação de Risco 


Consequência 


Decisão de projeto 


Figura 2.9. Passos de um estudo probabilístico. Modificado Haldar e Mahadevan (2000). 


Na avaliação de confiabilidade, os conceitos de estados limites são utilizados para 
ajudar na definição falha. Estes limites são considerados a fronteira entre os desempenhos 
desejados e não desejados para uma estrutura. Este limite é comumente representado 
matematicamente como uma função de desempenho ou função de estado limite. A teoria 
de probabilidade fornece bases precisas para medir a porcentagem de probabilidade em 
relação aos estados limites, assumindo um sistema em que suas variáveis sejam aleatórias 
e não mais determinísticas. Em casos mais simples, os estados limites são uma margem 
de segurança entre R (resistência ou capacidade) e Q (solicitação ou demanda). O estado 
limite da função g(.) é então definido em termos de R(.) e Q(.), que são funções de um vetor 
de variáveis aleatórias X, onde: 


IX) = R(X) - AX) (2.44) 


Com a probabilidade de falha, P,, expressa como: 


P,= PIg(X)<0] EM) 


Em que g(X)<o0 é a região de falha; g(X) = O é a superfície limite de falha; e g(X) >0 
é a região de segurança. 
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2.4.1 Conceitos de estatística 


Para quantificar as incertezas, é importante saber avaliar os riscos decorrentes 
da utilização de informações provenientes de observações, com base nos conceitos de 
teoria da matemática, a probabilidade, modelagem das incertezas em variáveis continua e 
discretas, determinação da distribuição e os parâmetros. Pode-se definir que as variáveis 
aleatórias são valores mensuráveis no espaço de probabilidade de experiências associadas 
com eventos. 

Ao analisar um evento, pode-se observar que não existem valores máximos e 
mínimos absolutos, pois os valores alteram-se de acordo com a quantidade dados. Para 
melhorar os problemas de valores máximos e mínimos, trabalha-se com o cálculo de uma 
média. Desta forma, utiliza-se um intervalo variável X, com n observações, onde a média é 
uma medida de tendência central dos dados, escrita como |, que pode ser calculada para 
n observações, conforme Equação (2.46). 


1 n 
EO == 4 (2.46) 
i=1 


Amédia (u,) pode ser expressa em termos de variação, desvio padrão e o coeficiente 
de variação (CV). A variação de X é denominada também como segundo momento central 
que é expresso: 


n 
1 
= > ==, P (2.47) 
i=1 


O desvio padrão expresso como o, fica nas mesmas unidades da média e é 
considerado a medida mais comum da dispersão estatística. 

O coeficiente de variação adimensional consiste em uma relação entre o desvio 
padrão e a média. 

CVO) = 8= E (2.48) 
x 

O diagrama de frequência, representado por um histograma, tem como objetivo 
marcar a forma da curva modelo que servirá para a determinação da Função de Densidade 
de Probabilidade (FDP) representada como f(x) como mostrasse na Figura 2.10. 

Pode-se também complementar esta informação, representado graficamente a 
informação em um histograma. A área abaixo da curva dependerá da largura dos intervalos 
e do número de pontos designados, pela probabilidade do evento acontecer entre O e 1. 
Matematicamente, isto é muito vantajoso porque a área total é igual a uma unidade. Assim, 
um histograma de unidade 1 é denominado como diagrama de frequência que é obtido por 
meio do diagrama de frequência ordinário dividido pela área. 
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Outra maneira de descrever a distribuição de probabilidade para as variáveis 
aleatórias discretas e contínuas é a Função de Distribuição Acumulada (FDC) ilustrada na 
Figura 2.11, F 4x). O FDC é definido para todos os valores de variáveis aleatórias X de -o 
ate oo que é igual à probabilidade de X ser menor ou igual a um valor x. Para uma variável 
aleatória contínua, F,(x) é calculada integrando a FDP para todos os valores de X menores 


ou iguais a X. 


FeO)= [5 fx (30) dx (2.49) 


Leto) 


Figura 2.11. Função de distribuição acumulada (FDC) 


O FDC das variáveis aleatórias discretas, pode ser matematicamente expresso 


assim: 


FX) =P(X <x) =Daxiex DX (Xi) (2.50) 
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Haldar e Mahadevan (2000) discutem em detalhe estes conceitos de estatística e 
ademais escreveram o livro de risco e avaliação de confiabilidade. 


2.4.2 Tipos de Distribuições Continuas 


Comumente, as variáveis aleatórias discretas e continuas são estudadas 
separadamente e, além disso, buscam fazer uma classificação do tipo de distribuições. A 
seleção da função de distribuição é muito importante para conhecer as características de 
probabilidade de um sistema estrutural. 


Distribuição normal ou Gaussiana 


Gauss (1800) estudou a distribuição normal ou Gaussiana, sendo esta uma das 
mais comuns e usadas em problemas estatísticos de engenharia, onde o FDP pode ser 
calculado como: 


mea = (rue? 
fe09 = x exp | à. (258) | -00 < x <+ 00 (2.51) 

Na Equação (2.51) e (2.52) tem-se que: yu, é a média e 0, o desvio padrão. O FDC 
pode ser calculado como se expressa na Equação 2.58. 


Fy(x) = IN : sete [5 (| dx (2.52) 


A distribuição normal é escrita como N (uy, 0) e tem características que facilitam 
o trabalho por suas propriedades, sendo aplicada para intervalos variáveis de -co até 00, 
Esta distribuição é simétrica com respeito ao valor médio. Também os valores modais são 
idênticos e podem ser designados diretamente dos dados. Nem sempre é possível estimar 
a probabilidade integrando a Equação (2.51) e nesse caso, é necessário a transformação 
das variáveis aleatórias X em variáveis normal padrão, com média zero e desvio padrão 
igual a 1, denotado como N(0,1). 


(2.53) 
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fel) 


[eo] 


Figura 2.12. Distribuição normal 


Após transformar as variáveis em normal padrão, a Equação (2.51) pode expressar 
o FDP de s como: 


fs(s) = == em (=:s?) -00 <s<+ 00 (2.54) 


Para o gráfico de FDC de S está expressa: 
1 


Fo(s) = [e (- 55?) ds (2.55) 


O FDC também é escrito como à(s), onde &(s) = F(S) Equação (2.55). A distribuição 
normal padrão é simétrica, como pode ser mostrado na equação a seguir: 


Ps) = 1.0 - ds) =P, (2.56) 


Distribuição lognormal. 

A distribuição lognormal foi proposta por Galton (1879). Em algumas situações, 
trabalha-se apenas com valores positivos devido à natureza do problema. Neste caso, é 
melhor utilizar distribuição lognormal que limita as variáveis aleatórias positivas e elimina 
qualquer possibilidade de valores negativos. Para distribuição lognormal, tem-se logaritmo 
natural com uma distribuição normal, ou seja, x é lognormal se y = Ln(x) fosse normal. 
Então, para uma distribuição lognormal o PDF está expresso da seguinte maneira: 


2 
fu) = —>exp|- +. (22) | g<x<+0 (2.57) 
oyV2nx 2 oy 
Onde: 
2 
oy = In (E) + | (2.58) 
X 
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=L 1 2 
co Gin iiS o da (2.59) 


Da Equação (2.57), e uy, e o, são dois parâmeiros da distribuição lognormal, 
denominados média e desvio padrão, respectivamente. As variáveis lognormal têm valores 
entre O <x< +º0, com FDP assimétrico, onde a média e os pontos modais esperados são 
diferentes da distribuição normal. 


0,25 


— 1=1,6=1 


sumo 4=2,0=1,5 


meses U3,0=] 


f(x) 


Figura 2.13. FDP lognormal 


Os 0,e y, são o desvio e a média equivalente e o FDC da distribuição lognormal é 
denotada pela Equação (2.60). 


(2.60) 


Fy(O) = cexp|-5 RL iea | 


El, 


Distribuição Uniforme. 


Nesta distribuição, a função de densidade é “plana” e contínua, portanto a 
probabilidade é uniforme em um intervalo fechado [A*, B*]. 
Para o intervalo [A*, B*], a variável uniforme contínua X dentro do intervalo está 


definida como: 


1 
fla, BN lp q4 <í<B (2.61) 
0 em outro caso 
A função de densidade da Figura 2.14 apresenta uma forma retangular com altura 


constante e base B* - A*. 
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A* Br 


Figura 2.14. FDP uniforme. 


A média e a variância são, respetivamente: 


A* + B* 
== (2.62) 
(B* — AM)? 
2 — 
qe (2.63) 


Distribuição Gamma. 


A distribuição gamma é uma função de densidade que tem muitas aplicações, como, 


a exemplo, a distribuição exponencial que é um caso particular da distribuição gamma. 
Esta distribuição é muito usada em sistemas elétricos com tempos de operação variáveis. 
A distribuição gamma deriva-se da função matemática gamma. 


A função gamma é definida por: 


T() = fi. xº-1e-* paraÃ>0 (2.64) 


Com a FDP da distribuição gama apresentada na Figura 2.15 e com os parâmetros 


A ey a variável aleatória X tem uma distribuição gamma se a função de densidade está 
denotada por: 


i 
f(x,4,7) =4y1T (4) 
0, em outro caso 
Onde A>0 ey >0. 


xe -XY x >0 (2.65) 


A média e a variância da distribuição gamma é definida por: 


v=N (2.66) 


2=Ay (2.67) 
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Figura 2.15. FDP gamma. 


Distribuição t. 

A distribuição normal trabalha com o conhecimento da média e o desvio padrão da 
população. É mais relevante saber o valor da média |, que o valor do desvio padrão o, 
porque com o desvio padrão pode ser assumido nos caso que tenha um amplo conhecimento 
do processo das amostras estudadas. Além disso, o valor o também pode ser estimado a 
partir do valor médio amostral. Para yu, tem se que: 


=: + (2.68) 


—“SNhn 


S, é análogo à o, na Equação (2.68). A distribuição t desvia-se de forma parecida à 


distribuição normal. Para amostras grandes com n > 30, a distribuição t é muito similar à 
distribuição normal. Para amostras com n< 30, é recomendável trabalhar com a distribuição 
t Dessa maneira, para uma distribuição t de uma amostra extraída de uma população 
normal, com v, graus de liberdade pode escrever-se como: 


po EBROND. A O (2.69) 


622 AWIA-D) 


Com uma distribuição normal padrão: 


n— DS 
po PD (2.70) 
q? 
Sendo V uma variável aleatória que segue uma distribuição X? chi quadrada com v, 
=n-1 graus de liberdade. 
É com £, Z, V, S? variáveis independentes, onde Z é uma variável aleatória normal 


padrão. A variável T pode ser expressa como: 
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Z 
T= 7 (2.71) 


E sua FDP é dada por: 


24 —(vy+1)/2 
TI(v.+1)/2] (1 t ) 


MO = To fr 


,—00 <t<oo (2.72) 


Vx 


A Figura 2.16 apresenta a distribuição t com v, graus de liberdade. 


t—a=—ta 0 ta 


Figura 2.16. FDP. 


Distribuição t-local com fator de escala 


A distribuição t-local com fator de escala é útil para modelar distribuições de dados 
com caudas mais pesadas que a distribuição normal. A distribuição t-local escala aproxima- 
se da distribuição normal como v, se aproxima do infinito, e para menores valores de v 
produzir caudas mais pesadas. 

A função de densidade probabilidade da distribuição t-local escala está expressa: 


ND W+D/2 
rico + 0/2) (vc + (EE) 
oTv(ve/2) Vx 


Onde T é a função gamma, | é um parâmetro local, o é um parâmetro local e v, é 


(2.73) 


um parâmetro de forma. 
A média da distribuição t-local escala é o parâmetro local y. Onde A média é definida 
apenas para valores de parâmetros forma v,> 1. Para os outros valores de v,, a média é 
indefinida. 
A variância da distribuição t-local escala está definida: 
Vx 
Ve — 2 


02=9?2 


(2.74) 
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Onde yu é um parâmetro local e v, é um parâmetro de forma. 


Valor extremo generalizado. 


Para a distribuição de valor extremo generalizado (VEG), a função de densidade 
probabilidade com os parâmetros |, o e parâmeiro de forma k, pode ser expressa como: 


1 1 
x “(1 (x — w)Y Fx (x— NR (2.75) 
y=[Clkoti0) = (5) exp (14%, - ) (1, - ) 


Para k, > O corresponde ao caso Tipo Il, para o fator k <O correspondente ao caso 
Tipo Ill e k, =0 caso Tipo |, onde a função de densidade está expressa: 


1 — = 
y = fCxl0,u,0) = (=) exp [esp (- ED) A Cd 


[0] 


(2.76) 


Sendo k, o parâmetro de forma que mede a dispersão. 

A distribuição de valor extremo generalizado a é utilizada para eventos nos quais 
se tem um valor más pequeno o um maior valor entre um conjunto de valores aleatórios 
independentes. 

O valor extremo generalizado mistura três distribuições simples, permitindo uma 
categoria continua de possíveis formas incluído três das distribuições mais simples. 

Pôr o teorema do valor extremo da distribuição GEV é a distribuição limite normalizada 
de máximos de um conjunto de variáveis aleatórias independentes e uniformemente 
distribuídos. A distribuição GEV utiliza-se como uma aproximação para modelar os máximos 
de eventos de variáveis aleatórias. 

Os três casos da distribuição de valor extremo generalizado fazem referência aos 
tipos |, Il, y Ill. As distribuições cujas caudas diminui exponencialmente como normal, 
levando ao tipo |. As distribuições cujas caudas com diminuição polinomial, tais como t 
de Student descrevem a tipo Il e as distribuições cujas caudas são finitos, tais como beta, 
levando ao tipo III. 

Tipos |, Il, y Ill fazem referência as distribuições de Gumbel, Frechet, e os tipos de 
Weibull. O Tipo | (Gumbel) e Tipo Ill casos de (Weibull). 


Revisão bibliográfica 


26 


ke <0, Tipo Ill 
o — = k =0, Tipo | 
k:>0, Tipo ll 


Figura 2.17. FDP Valor extremo generalizado. 


Distribuição Weibull. 


Weibull (1930) criou a distribuição Weibull, utilizada para sistemas complexos que 
dependem de vários componentes. Os componentes que constituem o sistema podem 
falhar por condições ambientais ou por outros parâmetros, que são difíceis de prever. 

Com os parâmetros À e y a variável aleatória contínua X apresenta uma distribuição 
de Weibull, se a função de densidade é definida por: 


Ayxr- te 
AW= ( >0 (2.77) 
654,7) 0, em outro caso Fá 


ParaA>0ey>0. 
A média e a variância para a distribuição weibull definem-se, respectivamente: 


1 
u=AYT (1 + 5) (2.78) 


o = an [r(145)-[r(1+5)]] 279 


A Figura 2.18 apresenta a FDP que descreve a distribuição Weibull. 
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Figura 2.18. FDP Weibull 


2.5 MÉTODOS DE CONFIABILIDADE ESTRUTURAL 


Choi, Grandhi e Canfield (2007) descrevem que existem muitos métodos para 
fazer análises de confiabilidade, citando-se: segunda ordem, simulação de Monte Carlo, 
Amostragem Hipercubo Latino e expansão estocástica. A falha estrutural e avaliação de 
confiabilidade são determinadas para os casos em que o estado limite da função não seja 
excedido. Estes métodos estão vinculados ao conhecimento da probabilidade de falha 
e, além disso, é importante conhecer as propriedades estatísticas, tipo de distribuição e 
os intervalos confidenciais de resposta das estruturas, que tem um papel fundamental na 
avaliação. 


2.5.1 Método FOSM 


A formulação original do método FOSM foi escrita por Cornell (1969). O Método 
de primeira ordem segundo momento simplifica o cálculo da probabilidade de falha. Esse 
método foi desenvolvido por meio da função de primeira ordem de expansão, com entrada 
e saída das médias e os desvios padrão. Nesse método, os momentos mais altos que 
descrevem a inclinação da distribuição são ignorados. 

No MVFOSM, a função de estado limite é representada como a primeira ordem de 
expansão em série de Taylor no ponto de valor médio, assumindo que as variáveis X são 
estatisticamente independentes. A aproximação do estado limite da função e a média são 


expressas como: 


HM = 949 + Vau) (X,- Us (2.80) 
Então, uy = (Ux, Hxo xo - - - Hx)", onde Vg(u) é o gradiente de g avaliado em ut. 
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(2.81) 


Ir) | 9Qtx) ao 
x, dx, "9x 


Vg(ux) = ( 
O valor da média avaliado no estado limite da função é expresso da seguinte forma: 
ug = Elg(ux)] = g(ux) (2.82) 


Devido a que a variância do estado limite e a variância do gradiente da função são: 


2 — 2 ne 
(gu = 0, vg] = O (2.83) 
o = aê — o? 

[Vou Tamo]  AvaqueTx]  Hx9vgq)] (2.84) 
= seê 
= 9Ivg(ux Tx] (2.85) 
= Vg(ux)” Po (2.86) 


O estado limite da função g(u,) tem uma variância que é expressa a seguir: 


2 As ad 2 
Taco * gw) + vamo” (x-mo] (2.87) 
O desvio padrão da função da aproximação do estado limite se dá por: 


E 
(É 2 
9 geo] = JIiBuo = [179 Gui)" oã = = [Er ( 29) * ) ol (2.88) 


Onde o índice de confiabilidade é: 


p= (2.89) 


Nos casos em que o estado limite da função seja não linear, é preciso fazer uma 


linearização do estado limite original da função ao ponto médio. Este método é conhecido 
como o método de valor médio MVFOSM e P é o índice de confiabilidade. Assim, para um 
espaço n-dimensional com variáveis independentes, a superfície de falha é hiperplana e 
pode ser expressa como uma função linear de falha pela Equação (2.90). 

n 


GM) = co + 5 CiXi (2.90) 
= 


O índice de confiabilidade para as n dimensões é encontrado com as Equações 
e (2.92). 
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Hg = Co + Cillxs + Collxa +": + Cnllyn (2.91) 


(2.92) 


O método MVFOSM simplifica os problemas complexos de confiabilidade com a 
utilização do índice de confiabilidade da Equação (2.89) que está relacionada aos valores 
de médias e desvios padrão. Ainda, esse método tem dois inconvenientes. O primeiro é que 
a linearização do estado limite da função ao valor médio pode levar a erros em funções com 
alta não-linearidade ou para coeficientes de variação grandes. 

O segundo inconveniente é que o método MVFOSM pode ser invariável em 
diferentes formulações matematicamente equivalentes do mesmo problema. Isto não é só 


um problema para formas não lineares de g (.), mas também para certas formas lineares. 


2.5.2 Índice de confiabilidade de Hasofer e Lind (HL) 


Hasofer e Lind (1974) propuseram um método que avalia g(X) em um ponto x* no 
espaço das variáveis aleatórias que se encontram na superfície do estado limite, como 
ilustra a Figura 2.19, para uma função linear g (R, Q) de duas variáveis aleatórias Re Q na 
forma da Equação (2.44). O ponto de superfície proposto para avaliação de g (R, Q), neste 
caso, tem as coordenadas (r*, q”) e é chamado ponto de projeto (PP). As coordenadas 
do PP são determinadas por um procedimento, onde as variáveis aleatórias sofrem uma 
transformação linear para um espaço de variáveis padronizadas e independentes Z. Um 
índice de confiabilidade, B, é interpretado como a menor distância da origem do espaço 
transformado até a superfície limite da função g(Z). No caso de funções g(Z), não lineares, 
como ilustra a Figura 2.19b, é necessário o uso de procedimentos de otimização matemática 
para achar a distância B através de algoritmos iterativos. 


Sra 


(a) Linear 
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(b)Não linear 


Figura 2.19. Funções de estado limite, índice de confiabilidade B e pontos de projeto para (a) função 
linear g (R, Q) e (b) função não linear g (Z) no espaço Z. Adaptado de Nowak e Collins (2013). 


Pode ser demonstrado que o índice £ se relaciona com a probabilidade de falha, P, 
através das seguintes relações: 


B=4-1(P)Ou 
P,= d(8) (2.93) 


Em que q é a função de distribuição normal padronizada acumulada. 
Como exemplo, considere a força R e a resistência Q, que são normalmente 
distribuídas e em que as médias y,, 1, € os desvios 0,, o, de Re Q, respetivamente. 


Então, as variáveis com normalização padrão são: 


R = tr q = Sto (2.94) 
OR IqQ 
Posteriormente, transforma-se a superfície do estado limite g(R, Q), com coordenadas 
originais do sistema (R, Q), dentro do estado limite de superfície nas coordenadas do 
sistema padrão normalizadas (À, 0). 


9 (R(R), Q(0)) =9(R,0)=Ror-Q0+(ur-uo)=0 (2.96) 


No sistema coordenado 5(R,0), a menor distância da origem até a superfície de 
falha 9(2,$) = 0, também conhecida como MPP, é o índice B = op+= (tecto) onde a 


tai A x ado DR oR2-0q? 
superfície de falha é a função não linear expressa por: 
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AX) = 94X, Xp XD) (2.96) 


Para uma distribuição normal padrão de uma variável u, os valores da média e o 
desvio padrão são O e 1. Logo, a transformação das variáveis utilizando q, e 0, é dada por:. 


Xi Uxi 
=== . 
i E (2.97) 


Região de segurança 
g>0 


Região de Falha 
Ed 


Figura 2.20. Geometria do index de seguridade. Modificado de Choi, Grandhi e Canfield (2007). 


A Figura 2.20 apresenta a superfície de falha do espaço X é definida no espaço 
U como g(U) = 0. O comprimento da origem no espaço U e o ponto g(U) = O é o desvio 


padrão do ponto médio no espaço X que corresponde ao ponto sobre g(X) = 0, onde B é o 
comprimento da superfície de falha (Equação 2.99). 


B(U) = (UTU)º = UM, VE g(U)=0 (2.98) 


u 
X> E 


U- Espaço 
X- Espaço 


Região de falha 
Região de 
segurança 


Região de 
segurança 


. Cy , Hxz) 


g(x)=0 54 


X1 


Figura 2.21. Superfície de falha do espaço X ao espaço U. Modificado de Choi, Grandhi e Canfield 
(2007). 
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Na Figura 2.21 ilustra-se o 8 que é comprimento mais curto entre a origem e à 
superfície de falha g(U) = 0. 


E 
Z 


B(U) = min(U'U) 
Ve gU)= 


(2.99) 


2.5.3 Método de iterativo de Hasofer e Lind (FOSM-HL) para distribuições 
normais 

Hasofer e Lind (1974) propuseram um método iterativo, denominado FOSM-HL que 
facilita os cálculos de um problema, para resolvê-lo por meio de iterações. No FOSM- 
HL, assume-se o estado limite da superfície com n-dimensões para variáveis aleatórias X, 
normalmente distribuídas e independentes. 

Inicialmente assume-se uma função de estado limite normalmente distribuída para a 
variável aleatória independente X como: 


9X) = 94X, Xp. XH) =0 (2.100) 
O estado limite pode ser determinado para funções lineares e não lineares. O estado 


limite da função descrito na Equação (2.100) pode ser transformado em normal padrão com 
a Equação (2.97) e ser expresso como: 


HU =9H0M, Vo SM Ho es SM tU STD =0 (2.101) 


NU = 960,4, + Vi SM FU e SM tUT)=0 (2.102) 


O índice B é o comprimento de MPP até a origem, em que a série de expansão de 
Taylor de primeira ordem vai de g(U) até o MPP U* e está expressa da seguinte forma: 


sa) = ur) + > da 2 (u ui — ui”) (2.103) 


Para a transformação das variáveis da Equação (2.97), designa a letra K para o 
número de iterações do algoritmo. 


dJ(U) 098 
du; [oba do 


(2.104) 


O índice de segurança é o menor comprimento da origem até a superfície de falha, 
onde O é o vetor normal da origem à superfície de falha g(U), que intercepta o ponto P*: 
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0,001. 


SU)-E, AD gui 


2 
n 
al Gui) 


(2.105) 


ôg(U") 
0x; 


dao ê) 
u; 
0. = s0,.= sie 
cos Ox; = Cos Oi Vo(UO] 
Te 49 
0X; XL 
gã=-———TI o — 7 (2.106) 
si, (EE ) É 
(PA Oxi 
As coordenadas para o ponto P* são calculadas assim: 
XxX — le A 
ui = pata o A OP* cos 0x; = B cos Oyi (2.107) 
Oxi 
As coordenadas do ponto P* em seu espaço original são: 
x*=4,+ BO, cos0,, (i=1,2,..,n) (2.108) 
Então PP* é um ponto na superfície do estado limite. 
GÃX, XX) =0 (2.109) 
Para saber se 8 converge, a desigualdade tem que ser satisfeita. 
= Bra Bi or (2.110) 


Bi 


Em que k é número de iterações e TOL é uma tolerância numérica especificada de 


As etapas do método (FOSM-HL) são as seguintes: 


1. 
2: 


definir a função de estado limite apropriado, da Equação (2.100); 


determinar o conjunto de médias para o ponto inicial de projeto x, = |, com 
i=1,2, ..., n. Em seguida, calcular os gradientes Vg(X,) da função de ásiado limite 
neste pon; onde os x,, representam os elementos do vetor X, e k representa 
o número de iterações; 


calcular o £ inicial usando o método do valor médio, B= H/0, e a direção do 
cosseno; 


calcular os novos pontos X, e U, das Equações (2.107) e (2.108) e os gradientes 
nestes novos pontos; 
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5. calcular o índice de segurança B usando a Equação (2.105) e a direção do 
cosseno com na Equação (2.106); 
6. refazer os passos (4) e (6) até que B convirja na Equação (2.110); e 


7. calcular as coordenadas dos pontos do projeto X,ou o mais provável ponto de 
falha (MPP), X*. 


2.5.4 Método Hasofer Lind-Rackwitz Fiessler (HL-RF) 


Para este método, as variáveis aleatórias X são consideradas normalmente 
distribuídas. Em casos não gaussianos, mesmo quando a função g limite de estado (X) 
é linear para o cálculo da probabilidade, precisa-se fazer uma transformação. Existem 
diversos métodos disponíveis para a realização de transformações, tal como e Rackwitz. 
Um simples, transformação de aproximação chamada a distribuição normal equivalente. 

Quando as variáveis são independentes entre si, onde q[. ] é a inversa de dI. ], a 
transformação é: 


U,= 4º [FX] (2111) 


Para obter uma distribuição normal equivalente, utiliza-se na série de expansão de 
Taylor para a transformação de MPP até X*, não levando em consideração os termos não 


lineares. 
ô 
u; = 4 "HFu(x)] + 34, 14 Fai Ce) D ir Qi =Xi) Ra 
Onde: 
FS 2113 
ER [Fxi(%i)] db HPa(x)]) | 


Substituindo o termo 4 EaD] da Equação (2.113) na (2.112), tem-se: 
MT li o 4 [El] (d Fi )D/ fel] (2.114) 
PCP Fxi(xi DD / Pei lxi!) 


epa dust) (2.115) 


Onde y,, e 0,, são, respectivamente, a média e o desvio padrão equivalentes ou 
aproximados da distribuição normal e F, (x,) éa FDC e f,(x,) é a FDP. 


— POTE] 


é Tas 2.116 
GT Fixo) Ea 
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Hei = x" — 0 Falow (2.117) 
Para obter distribuições normais equivalentes, trabalha-se com a FDC e com FDP 
da função original não gaussiana com variáveis aleatórias. Posteriormente, transformam- 
se as variáveis por meio de aproximação em variáveis aleatórias equivalentes normais ao 
MPP. Com x' sendo uma variável aleatória, com distribuição equivalente, normalmente 
distribuída. 


FAX) = FoolXT) (2.118) 
FaQi) = 9 (2.119) 
Hei =" — 4 Ec (2.120) 


Os passos do método Hasofer Lind — Rackwitz Fiessler (HL-RF) 

A presente pesquisa propõe a implementação da metodologia de Hasofer e Lind 
(1974) melhorando-as com o método de Rackwitz e Fiessler (1978) para realizar o 
tratamento de variáveis não normais. O algoritmo foi implementado com o processo iterativo 
de Newton-Raphson. O organograma Figura 2.22 demostra a organização do processo 


implementado. 


Inicialização 
Definir 900) = g(L,x2,...xn)) = 0 
i=1 
Conjunto de tolerâncias de convergência 4, » (e.g. 1073) 
Escolher um ponto de partida x, (e.g. x4= Hy) 
u, = u(x,) 
Parâmetro de escala Go(e.g. Go= Hx) 


x; = x(u;) 


Jul a) na ] 

* G(u;) = g(x;) Atualização 
* VG(uj)= Ve(xi) Jxu G(u) |aT 
* é =-V6(u)/VG(uo)| E tes [s+ E) e 


* Algoritmo Newton-Raphson modificado e 
melhorado para uma melhor convergência. 


Verificação 
verificar convergência 
* |G(u)/Grilses 

n lu; = Quê || < E2 


: Notas: : 
: xy=x(u,e Ju SÃO à transformação : 
; e O respectivo Jacobiano entre a | 
: distribuição de x e distribuição : 
; normal padrão não correlacionado ; 
: onde u-N(0,1) : 


Sim 


Finalização 

"“ âà=á; 
su =u; 

“ B=âu' 

“ x'=x(u”) 

- Pr=0(-B) 


Figura 2.22. Procedimentos para o cálculo do índice de confiabilidade B para funções de estado limite 
lineares e não lineares 
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2.6 MÉTODO DE MONTE CARLO (MMC) 


Segundo Haldar e Machadavam (2000), o método de simulação de Monte Carlo é 
composto pelos seguintes passos: 


* | fazer a formulação do problema e definir as variáveis aleatórias; 

* quantificação das características probabilísticas das variáveis aleatórias; 
* | geração dos números aleatórios; e 

* | experimentação numérica. 


A formulação do problema identifica o evento a ser estudado e as variáveis que 
estão relacionadas. 

Em seguida, para fazer a quantificação das características probabilísticas das 
variáveis aleatórias, define-se inicialmente o possível comportamento que apresenta uma 
população e que está definido pela função de densidade ou distribuição da probabilidade 
FDP. 


b 
Pla <x<b] = | f(x). dx (2.121) 


Em que F é a função de distribuição acumulada de X: 


Púdi= [ “ Fudu (2.122) 


Para variáveis que apresentam um comportamento normal Gaussiano: 


FDP = O 


1 
n2==). E da (2.123) 


A função de distribuição normal padrão é expressa: 


u2 


FDP = Poa(x) = -Z du (2.124) 


1 X 
—=).e 
ov21 Jo 

É possível relacionar todas as variáveis aleatórias normais com a distribuição normal 
padrão, se X- N(u, &º). Então: 


PR É (2.125) 
[0] 


Este processo denomina-se padronização X- N (yu, &) > Z- N(0,1) . Se Z é uma 


distribuição normal padrão com y e o sua inversa, então é expressa por: 


X=0Z+y (2.126) 


A simulação de MMC à FDC também possui um papel importante do processo, pois 
é utilizada nos cálculos para a geração de números aleatórios. 
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O processo de geração de números aleatórios, inicialmente, identifica as 
características probabilísticas das variáveis aleatórias, como as médias, os desvios padrão 
e o tipo de distribuição. MMC baseia-se em gerar números aleatórios de acordo como 
a distribuição especificada. A geração de números aleatórios é realizada para números 
aleatórios distribuídos uniformemente entre O e 1. Em seguida à geração dos números 
aleatórios uniformes u,, esses valores são transformados em uma distribuição normal 
padrão, ou em qualquer outro tipo de distribuição seguido o mapeamento descrito na Figura 
2.23 que é baseado no teorema da integral equivalente entre as FDC. O MMC se utiliza 
deste terorema para a geração dos números aleatórias necessária à simulação. 


Fu(u) Fx(x) 


á fu(u) fx(x) 


Figura 2.23. Mapeamento para simulação de números aleatórios (modificado de Haldar et al., 2000) 


A variável aleatória de FDC é igualada ao número aleatória gerado em ui. F(X,) = 
u 


xi? 


com a seguinte equação: 


X = Fa (u,) (2.127) 


A experimentação numérica é realizada para N simulações, até não obter muita 


diferença entre as respostas do sistemas para cada N simulação. 
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2.7 TRABALHOS DE CONFIABILIDADE E MECÂNICA DA FRATURA 


Alguns trabalhos no Brasil relacionados com a pesquisa são os artigos (Diniz et 
al., 2011; Diniz et al., 2013) onde estudam-se os benefícios de durabilidade em relação 
ao corrosão do aço e a falha dúctil para o concreto armado e a falha frágil para o concreto 
armado com polímeros de fibras de carbono. Diniz et al., 2013 considera as incertezas 
das variáveis do projeto de vigas estimadas com a norma ACI-440 (2006), e faze uma 
avaliação de probabilidade de falha com respeito à resistência à flexão última por meio 
do método de Monte Carlo. Os procedimentos apresentados no artigo podem ser usados 
para o tratamento e recomendações de vigas de concreto armado com polímeros de fibras 
de carbono. Diniz et al., 2011 estabelece comparações entre as vigas intactas e as vigas 
recuperadas, onde os deslocamentos das vigas são avaliados com um modelo não linear 
para o concreto armado utilizando o método dos elementos finitos e simulações de Monte 
Carlo para conhecer a estatística do processo de carga ruptura. 

O artigo de Beck et al., 2008 estuda tubos de concreto armado com abertura de 
fissuras, projetado pela norma NBR 6168:2003. Nesse trabalho, faz-se uma avaliação das 
formulações respeitavas para o projeto de elementos de concreto armado por meio da 
confiabilidade utilizando o método FOSM e com simulações de Monte Carlo. Os resultados 
revela a influência das variáveis de projeto que tem maior influência na probabilidade de 
falha, relacionadas com o problema e com os índices de confiabilidade. 

Outro artigo do mesmo autor Beck et al., 2004 desenvolve uma pesquisa para O 
problema de falha por sobrecarga de um elemento estrutural, submetido a um carregamento 
aleatório para o qual analisa-se a propagação de trincas. O problema é formulado com a 
mistura de modelos de falha e modelos de propagação de trincas. A solução do problema 
aborda a sobrecarga e os modelos críticos de falha por meio do cálculo de crescimento 
da trinca evitando controles destrutivos. Além de isso o estudo de a sensibilidade de 
participação das variáveis K fator de intensidade de tensão, a, comprimento crítico de 
trinca, o, entre outras, ajuda a determinar o método de solução apropriado para problemas 
específicos. 
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METODOLOGIA 


Neste capitulo, é descrita a metodologia utilizada para a avaliação de confiabilidade 
realizada para a flexão de vigas em três pontos (FTP) e elaboradas com concreto 
convencional (CONC) e com concreto reciclado (FRAP). 


Figura 3.1. Metodologia geral 


3.1 PROGRAMA EXPERIMENTAL 


Ensaios seguindo as normas RILEM (1990a, b) para vigas submetidas a FTP de 
diversos tamanhos foram testadas para o cálculo das propriedades de fratura (G, C, 
CTOD,) e módulo de elasticidade (E). Os ensaios foram realizados no Laboratório de 
Materiais Cimentícios da Universidade de Illinois em Urbana Champaign (UIUC), nos 
Estados Unidos, para o CONC. A Tabela 3.1 mostra as dimensões e os dados medidos das 
vigas ensaiadas para os ensaios de EEB, MFLE e MDP para CONC. 

É importante observar que as mesmas vigas foram usadas em cada metodologia 
de ensaio. Assim, o EEB, o MFLE e o MDP determinaram três parâmetros estatísticos de 
média aritmética (1), desvio padrão (0) e coeficiente de variação (CV) das propriedades dos 
materiais (G,, C, CTOD,). Os resultados do cálculo desses parâmetros são apresentados 
na seção 2.3, em que P,, como se ilustrou na seção 2.3.1, faz referência ao carregamento 
máximo medido no ensaio. 

Ademais, é importante notar que P é um parâmetro do ensaio EEB, também usado 
para o ensaio de MFLE e que C e C são parâmetros dados do ensaio de MDP. As variáveis 
de L, Det estão relacionadas com a geometria da viga, Sé o vão e Wé o peso da viga. 
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Tabela 3.1. Dimensões e dados medidos das vigas para os ensaios de EEB, MFLE e MDP para CONC. 


ga O Emi) a no Ei (mi Es) ara E 
la 530,7 11050 100060 2550 8060 880 67  53E06 88-06 
lb 527,3 11050 10060 26060 790 950 69  60E06 88-06 
lc 529,2 11040 10040 2580 790 885 62  53E06 68,1E-06 
la | 1923 7000 6010 1520 790 580 43  72E06 1,4E-05 
lb [1937 7000 6020 1530 790 530 51  65E06 12-05 
lc | 1898 7000 6010 1540 790 550 47 66E06 1,1E-05 
lla | 40,6 3490 2520 630 8060 225 25  64E06 1,1E-05 
mb [4140 3475 2490 635 8060 220 26  70E06 1,1E-05 
lc | 414,8 3500 2490 640 795 20 22 - : 

ld [41,9 3475 2480 655 795 25 24 80E06 12-05 


O concreto reciclado tem uma utilização muito frequente devido os benefícios 
ambientais e econômicos e é composto por materiais reciclados proveniente de subprodutos 
de concreto. Alguns materiais utilizam o concreto reciclado em sua composição, tais 
como o cimento pozolânico e escoria de alto forno de terra granulada. Aditivos minerais, 
como cinzas volantes, também conhecidas como agregados reciclados, fracionados ou 
recuperados de pavimento asfáltico são utilizados como agregado para confecção de 
concreto reciclado (RCA). Entre os principais utilizados estão o RAP do inglês recycled 
concrete aggregate e o FRAP do inglês fractionated reclaimed asphalt pavement. 

O concreto com adição de agregado reciclado denominado FRAP é elaborado da 
mistura de cimento com agregado obtido de pavimento asfáltico fraturado e reciclado, e 
tem propriedades de capacidade de carga de flexão de losas de concreto especificamente 
muito parecidas com a do concreto convencional CONC como observa-se nos resultados da 
pesquisa realizada por Brant et.al., 2013. Devido a isto e aos custos de produção menores 
comparados ao concreto convencional, a aplicação de FRAP e RAP em pavimentos de 
concreto de rodovias e aeroportos tem se difundido bastante. A Figura 3.2 apresenta as 
superfícies de um corpo de prova de FRAP. 
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Figura 3.2. Concreto Reciclado (FRAP). Brand et al (2013) 


Os mesmos ensaios realizados para o CONC, descritos anteriormente, foram 


realizados para vigas submetidas a FTP de diversos tamanhos com o FRAP baseados em 


ensaios que foram realizados no Laboratório de Materiais Cimentícios da Universidade 


de Illinois em Urbana Champaign (UIUC), nos Estados Unidos. A Tabela 3.2 mostra as 


dimensões e os dados medidos das vigas ensaiadas para FRAP. 


Tabela 3.2. Dimensões e dados medidos das vigas para os ensaios de EEB, MFLE e MDP para FRAP. 


viga l Emi) Ra Re Ei mi ad ind Ro 
la 523,9 11070 100060 2590 805 885 68  68E06 1,3E-05 
lb 540,5 11150 100060 2570 790 860 58 73E06 12-05 
lc 529,2 11090 10010 2600 7950 9040 57 6906 1,2E-05 
lb [1874 7,000 600,0 1520 790 530 36 86E06 1,3E-05 
lc | 1879 6990 60140 1510 8060 5060 37 68E06 1,1E-05 
lla | 40,0 3480 2490 630 8060 21,0 19 78E06 13-05 
mb [414 3490 2520 650 8060 190 22 84E06 1,6E-05 
lc | 398 3470 2500 640 8060 190 24 68E06 1,5E-05 


3.2 VARIÁVEIS DE ANÁLISE E FUNÇÕES DE DESEMPENHO 


A análise risco é realizada nessa pesquisa com a predição do índice de confiabilidade 


para a falha estrutural de uma viga FTP, como ilustrada na Figura 3.3, e é submetida ao 


carregamento no meio do vão (P). 
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Figura 3.3. Viga submetida a flexão em três pontos. 


A função de desempenho desta configuração é dada por: 


tD GE 


1 

z 

Gi(P, a9,E,G,,C) = = 
n 


——""1——— |) —P =pHEB-p 3.1 
apena) == em 
Em que, para a viga FTP, g(a,) está definido pelas Equações (2.10) e (2.20), C, pela 
Equação (2.24) e PEÉB é termo o proveniente da predição de resistência do modelo de 
EEB, definido na Equação (2.25) para EEB e que é a mesma função de desempenho para 


o MDP e, ainda, em que, daqui em diante, será referido como PFEB. 
1 
tD/ GE 12 Mie 
Go(Bo.00.E,67,0) = (a (a0)) — E = Pr -—B (3.2) 


Para a segunda função de desempenho, o primeiro termo do lado direito da equação 
(PMFLEY é proveniente da predição de resistência do ensaio de MFLE definido na Equação 
(2.26) para MFLE e que será referido como PMFLE 

O segundo termo das Equações (3.1) e (3.2) é o termo de solicitação (P) que será 
uma variável aleatória com média aritmética para cada «, e D, adotados como: 

1 
Hedao = FME so (3.3) 

Em que F, é um fator considerado aqui de F, =1,1; 1,5 e 2,0 e para CONC e de F. 
=1,1 para FRAP. 

Com G( )eG,( ) definidas pelas Equações (3.1) e (3.2), as propriedades estatísticas 
para o carregamento aplicado e os parâmetros geométricos são descritos na Tabela 3.3 
e na Tabela 3.4, respectivamente. Com os dados dos ensaios de laboratório e a teoria 
descrita nas seções 2.3.1. e 2.3.2, pode-se calcular as propriedade de fratura do material 
das vigas seguindo cada metodologia de teste (EEB, MFLE e MDP). As tabelas mostram 
os parâmetros estatísticos de média aritmética (4), coeficiente de variação (CV) e tipo de 
distribuição analisados. 

As análises de confiabilidade estudam dois casos de tipos de distribuição (normal 
e lognormal) para as variáveis aleatórias consideradas, como indicado da Tabela 3.3 à 
Tabela 3.6, para o Caso | e Caso Il para CONC e o Caso Il para FRAP. Estas distribuições 
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são inicialmente consideradas por serem as mais utilizadas na engenharia, especialmente 


quando há falta de testes extensivos com muitas amostras para determinar a respectiva 


distribuição. 


Tabela 3.3. Parâmetros de carregamento para as simulações de confiabilidade 


Variável u CV Caso | Caso Il 
Variável de acordo 
PAN) Com Equação (3.3) 0,10 Normal Lognormal 
Tabela 3.4. Parâmetros geométricos para as simulações de confiabilidade 

Variável u Cv Caso | Caso Il 
S (m) (0,25; 0,60; 1,00; 4,00) - Determinístico Determinístico 
t(m) 0,08 x Determinístico Determinístico 
D (m) (0,063; 0,15; 0,25; 1,00) - Determinístico Determinístico 
a, (m) Variável de modo que a,=+0,10; 0,10 Normal Lognormal 


0,25; 0,40) 


De acordo com os ensaios de laboratório considerando os dois procedimentos 


estudados EEB e o MDP, a Tabela 3.5 e 3.6 apresentam, respectivamente, as características 


estatísticas das propriedades elásticas (E) e de fratura (G,e C) dos materiais determinados 
com cada ensaio. 


Tabela 3.5. Parâmetros do material, determinados pelo EEB, para as simulações de confiabilidade 


Material CONC FRAP Caso | Caso Il 
Variável V Cv V CV 

G, (Nim) 70,29 0,17 73,85 0,33 Normal Lognormal 
C (mm) 37,10 0,01 50,32 0,016 Normal Lognormal 
E (Nm?) 34,82E+09 0,07 25,83+09 0,11 Normal Lognormal 


Tabela 3.6. Parâmetros do material, determinados pelo MDP, para as simulações de confiabilidade. 


Material 


CONC FRAP Caso | Caso Il 
Variável V Cv V CV 
G, (Nim) 61,74 0,16 56,55 0,33 Normal Lognormal 
C (mm) 28,92 0,19 31,76 0,38 Normal Lognormal 
E (Nim?) 34,82E+09 0,07 27,04E+09 0,06 Normal Lognormal 
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ANÁLISE DE CONFIABILIDADE EM VIGAS FTP 


Neste capitulo apresentam-se os resultados das análises de confiabilidade 
realizados na seção 4.1. para o concreto convencional CONC e na seção 4.2. para o 
concreto reciclado FRAP. 


4.1 CONCRETO CONVENCIONAL (CONC) 


Os resultados de confiabilidade e a predição determinística de resistência para 
diferentes tamanhos de viga FTP e de trinca, variando os a para material CONC avaliado 
com o algoritmo (HL-RF). 


4.1.1 Predição determinística do efeito de escala da resistência em vigas FTP 


Uma análise determinística considera a predição da resistência nominal (0,) com 
relação à geometria da viga para os três modelos estudados: EEB, MFLE e MDP. A Figura 
4.1 apresenta a dependência de o, com o tamanho das vigas (D) e o tamanho relativo das 
trincas (a, = c,/D) para o fator de segurança F. = 1,1. Os gráficos foram gerados de acordo 
com a análise determinística usual, em que são utilizadas as médias (u) dos parâmetros 
das Tabela 3.1, 3.3, 3.4, 3.5 e 3.6 para CONC. 

Como esperado, o efeito de escala no concreto é descrito pelo EEB por meio da 
diminuição da resistência nominal com o aumento do tamanho da viga, como ilustram as 
tendências das curvas da Figura 4.1 para os itens (a) e (b). Como esses dois modelos 
utilizam a mesma equação de predição para o,, as curvas de predição são bem similares e 
as pequenas diferenças podem ser atribuídas às diferenças dos parâmetros G,e C que são 
maiores quando calculados segundo a metodologia EEB como encontra-se especificado 
na Tabela 3.5. e na Tabela 3.6. É importante observar que, para a> 0 na menor viga (D 
=0,063m), os modelos EEB e MDP apresentam falhas de predição devido ao fato de que 
o tamanho equivalente da zona de processo ser a propriedade C, que nestes casos, tem 
o valor de C (28 mm/37 mm), e aproxima-se significativamente da espessura da viga (63 
mm). Este fato torna limitadas as análises nestes casos devido as hipóteses básicas dos 
modelos EEB e MDP (Bazant; Planas, 1997; Shah et al., 1994). 

As predições da MFLE, na Figura 4.1 para o item c), são bem conservadoras quando 
comparadas às preditas pelo EEB e MDP. Isto ocorre especialmente para pequenos valores 
de a, onde o, cresce assintoticamente na Equação (2.8), que também pode ser definida em 
termos da Equação (2.26). 


Análise de confiabilidade em vigas FTP 


45 


o (N/m?) 


o (N/m?) 


3,5E+06 


— (11=(),1 — (12=(),25 — (13=(),4. 


3,0E+06 


2,5E+06 


2,0E+06 


1,5E+06 + 


1,0E+06 


5,0E+05 


0,0E+00 7 7 
0,6 0,8 1 


3,5E+06 + — (110,1 —(2=0,25 — (13=(),4 


3,0E+06 + 


2,5E+06 + 


2,0E+06 + 


1,5E+06 + 


1,0E+06 + 


1 


5,0E+05 


0,0E+00 r 7 


fo 0,2 0,4 0,6 0,8 1 
D (m) 


(b) MDP 


Análise de confiabilidade em vigas FTP 46 


8,0E+06 


7,0E+06 


6,0E+06 


5,0E+06 


4,0E+06 


e (N/m?) 


3,0E+06 


2,0E+06 


1,0E+06 


0,0E+00 


— (11 =(),1 — (12 =(),2 5 — (13=(),4 
« 
T T T T 1 
0,2 0,4 0,6 0,8 1 
D (m) 
(0) MFLE 


Figura 4.1. Predição deterministica da resistência o, e efeito de escala para a viga de FTP (CONC): (a) 


EEB; (b) MDP e (c) MFLE. 


A Figura 4.2 é referente ao gráfico clássico do efeito de escala que descreve a curva 


da resistência nominal com o aumento de tamanho em escala logarítmica como ilustrados 
na Seção 2.3.1. Utilizando a Equação (2.11) e (2.13) pode-se observar que as curvas são 


as esperadas, conforme ilustrado na Figura 2.5. 


10 


+) 


log(o,/Bf, 


a2=0,1 —— 43=0,25 —— q4=0,4 


1 10 
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a2=0,1 a3=0,25 —— q4=0,4 


log(o,/Bf 


1 10 
log (D/Do) 


(b) MDP 


Figura 4.2. Predição deterministica da resistência o, e efeito de escala para a viga de FTP em escala 
logaritmica para (CONC): (a) EEB; (b) MDP. 


4.1.2 Predição do índice de confiabilidade Pp e efeito de escala da resistência 
em vigas FTP 

Para a análise de confiabilidade, os índices B foram determinados com o algoritmo 
descrito na Figura 2.22 para o problema definido na Equação (3.1) e (3.2), com as condições 
das Equação (3.3). A Figura 4.3 apresenta os resultados de B para diferentes geometrias 
de viga e tamanhos relativos de trincas (a= a,/D) para F, = 1,1. Como indicado nessa figura, 
as análises foram realizadas para as propriedades extraídas dos diferentes modelos EEB, 
MDP e MFLE para o caso |, que definem as variáveis aleatórias seguindo uma distribuição 
normal (ver Tabela 3.1, 3.3, 3.4, 3.5 e 3.6). 

É importante notar que a função de desempenho G = P,- P faz relação entre os 
modelos de EEB, MFLE e MDP, trabalhando P, para cada modelo em particular (EEB, 
MFLE e MDP e, no casso de P,, trabalha-se com o modelo de EEB. Fazer estes tipos 
de comparações incrementado o valor de estado limite da função e os valores do índice 
de confiabilidade, mas, em contrapartida, é mais difícil a linearização e convergência, 
especialmente para valores pequenos D = 0,063 m e de a > 0,0 como fosse mencionado 
anteriormente. Portanto, compara-se EEB com todos os modelos, com exceção de MFLE 
para a> 0,0. Esse é comparado com o próprio para conseguir a convergência. 
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Figura 4.3 Predição do índice de confiabilidade B para material CONC, com F,= 1,1 e caso |: (a) EEB; 
(b) MDP e (c) MFLE 


Similarmente, a Figura 4.4 apresenta as predições de confiabilidade para o caso 
Il em que as variáveis aleatórias envolvidas são assumidas lognormais, com F. = 1,1. As 
mesmas tendências observadas para o caso | podem ser aplicadas quando as variáveis 
são lognormais. Somente as predições da MFLE mostraram um ligeiro aumento dos valores 
de B. A consideração da lognormalidade das distribuições de variáveis aleatórias é muito 
importante uma vez que variáveis aleatórias lognormais não assumem valores negativos. 
Isto se torna muito conveniente na consideração de variáveis de geometria e resistência 


comumente encontrados em problemas de engenharia. 
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Figura 4.4. Predição do índice de confiabilidade B para material CONC, com F, = 1,1 e caso Il: (a) 
EEB; (b) MDP e (c) MFLE. 


Os índices 8 forem muito maiores para MFLE em comparação com os modelos de 
MFNLE de EEB e MDP, devido a que MFLE não considera a incidência do parâmetro de 
comprimento de trinca equivalente e além disso não consegue convergir para pequenos 
valores de D e de a. É importante observar os altos valores para os índices de confiabilidade 
preditos pela MFLE, onde, segundo este modelo, a viga teria uma falsa probabilidade de 
falha bem pequena o que é são predições não realistas. Isto alerta o fato da limitação 
dos usos dos princípios da mecânica da fratura linear para materiais quasi-frágeis como 
o concreto, onde a existencia de uma zona de processo, aqui representada por C, altera 
significante as predições de falha dos materiais. 
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Comparando os índices calculados para o EEB, MDP e MFLE, nota-se que devido 
ao efeito de escala do problema, os índices 8 decrescem com o aumento da espessura da 
viga D. Isto evidencia o efeito de escala também nos índices de confiabilidade, onde vigas 
mais espessas tem menor confiabilidade e, consequentemente, maior será a probabilidade 
de falha comparando-se com vigas de menor altura. Note que o EEB apresenta o menor 
efeito de escala para B, como ilustram as curvas quase constantes em D, especialmente 
para os menores a. Isto pode ser explicado pelo pequeno valor do CV=0,016 de C (Tabela 
3.5) determinado pelos ensaios de EEB, que é significantemente menor que a variação 
de C obtida pelos ensaios de MDP. A análise da influência de C será comprovada logo em 
seguida. 

Conforme esperado, os valores crescentes do tamanho da trinca relativo à altura, 
dados pela variável a, diminuem a confiabilidade, aumentando, assim, a probabilidade de 
falha. Seguindo a mesma tendência dos gráficos de predição de resistência, percebe-se 
que os valores de B determinados pela MFLE são maiores (predições não conservadoras 
do risco de falha) devido à predição não conservadora de o, dado por esse modelo. Isto 
evidencia as limitações da MFLE na predição quantitativa de falha em estruturas de 
concreto. 

Comparando os modelos EEB e MDP, nota-se que o MDP apresenta predições de 
confiabilidade maiores (mais conservadoras) para todos os tamanhos de viga De tamanhos 
relativos de trinca a. Por exemplo, a curva a =0,1 apresenta 0,75 < B < 0,30 na variação em 
D. Por outro lado, o EEB apresenta variação 0,60 < B < 0,55 para o mesmo a. 

Para os fatores de segurança F, = 1,5 e F,= 2,0, apresentados no Anexo A e B. 
Um resumo dos índices os índices B para os diferentes F, são apresentados na Figura 
4.3 para um a =0,25 para o material CONC para distribuição normal (caso |). Note que as 
curvas têm um comportamento similar ao apresentado para F,= 1,1, mas com índices de 


confiabilidade maiores e com curvas tendendo a linearidade. 
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Figura 4.5. Índices de confiabilidade CONC para Caso | com a=0,25 e FS = 1,1;1,5 e 2,0 para: (a) 
EEB; (b) MDP e (c) MFLE 


A Figura 4.6, Figura 4.7 e Figura 4.8 apresentam os gráficos em 3D para o material 
CONC e caso Il, analisando o comportamento de os indices 8 em função de ae D, e 
estabelecendo comparações entre os modelos EEB com o MDP para os fatores F, =1,1; 
1,5e2,0. 


D(m) io 


a 


Figura 4.6. Índices de confiabilidade 8 com Fs = 1,1 para os modelos EEB e MDP, material CONC 
(caso II) 
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D (m) o 


Figura 4.7. Índices de confiabilidade B com Fç=1,5 para os modelos EEB e MDP, material CONC 
(Caso Il) 


D (m) So 


Figura 4.8. Índices de confiabilidade B com F, = 2,0 para os modelos EEB e MDP, material CONC 
(Caso Il). 


As superfícies de 8 para os modelos de EEB e MDP e os fatores de F, =1,1, 1,5 
e 2,0 ilustram a incidência do efeito de escala, onde os valores de 8 aumentam com o 
incremento da altura de viga D, e diminuem, significantemente, com o aumento de a. Além 
disso, as superficies geradas dos indices 8 para EEB são inicialmente menores que as do 
MDP com F, =1,1, e com o aumento de F. para 1,5 e 2,0 as superfícies para o modelo EEB 
apresentam maiores índices de confiabilidade que MDP. 
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4.1.3 Verificação do algoritmo implementado por (FORM HL-RF) com o Método 
de Monte Carlo (MMC). 

O algoritmo de confiabilidade (HL-RF) foi comprovado, por meio de simulações do 
Método de Monte Carlo (MMC), para os resultados do material CONC com o fator F. =2,0 
seguindo as mesmas equações de desempenho de estado limite com G, ( )e G,( ) para 
os modelos EEB, MDP. A escolha do fator 2,0 é devido a este gerar maiores valores de 
probabilidade de falha, sendo o caso mais crítico para as predições do algoritmo FORM 
HL-RF. 

Foram usados os mesmos parâmetros estatísticos e determinísticos do material 
estabelecidos no capitulo 3 para o caso | e utilizando N,,,=9E+7 simulações do MMC. 
Utilizando as equações de desempenho em conjunto com o MMC, temos como probabilidade 
de falha: 

NºG()<0 


P 
f Numc 


(4.1) 
Em que P, é a probabilidade de falha obtida de número de vezes que função de 


estado limite seja menor que O, pela divisão do N,, 


mo dU€ é O número de simulações de 


MMC. Então conhecendo P, é possível deduzir o índice Luc, Que é o índice de confiabilidade 
calculado pelo MMC, da Equação (4.3) com a Equação (4.2): 


P= PIB vel (4.2) 


Bymo= " PLPÍ (4.3) 


Em que q é a FDC da distribuição normal padronizada. 


— 1Bmmc-Brorull 
Bmmc 


pela Equação (4.3) e os B,,,y calculado pelo algoritmo implementado e, cujos resultados 


A Figura 4.9 apresenta a diferença relativa Ber , entre os B,mo calculados 
foram apresentados nas seções anteriores. Note que, para os dois modelos, as diferenças 


relativas, 8. , encontram-se abaixo de 0,04 (4%) em geral. Mostrando a acurácio do 


er? 
algoritmo FORM implementado. Somente para o MDP e com a > 0,4, cujos coeficientes de 
variação das variáveis de entrada são bem maiores, os valores de f,, tendem a aumentar, 


mas ainda abaixo da diferença relativa de 0,1. 
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Figura 4.9 Erro relativo B, entre os índices de confiabilidade calculados entre o algoritmo FORM e o 
MMC para F, = 2,0, material CONC, Caso | e para os modelos a) EEB e b) MDP. 
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4.1.4 Índice de sensibilidade á, das variáveis na determinação da confiabilidade 
para CONC 

De acordo com o algoritmo implementado e ilustrado na Figura 2.22, a quantidade 
&; é calculada por: comparações 


Q;= -VG(u*IVG(u*) II (4.4) 


De acordo com a Equação (4.4), os &; são os cossenos diretores indicando a 
direção do ponto mais provável de falha (u*), que está a uma distância 8 da origem no 
espaço normal padronizado. Assim, a magnitude de cada componente &;, sendo io índice 
relativo a cada variável, é um indicador direto da importância daquela variável para o índice 
de confiabilidade. 

Os valores de à; para o caso | e F, =1,1 são apresentados da Figura 4.10 à 
Figura 4.12, para os modelos EEB, MDP e MFLE, respectivamente. As mesmas figuras 
para F, = 1,5 e 2,0 são apresentadas no Anexo B. Os gráficos consideram diferentes a 
(= 0;0,1; 0,25; 0,40) e D (0,063; 0,15; 0,25; 1,0 m) nas análises de confiabilidade dos três 
modelos estudados (EEB, MDP e MFLE). Nota-se que, como explicado anteriormente, a=>0 
não pôde ser simulado para a MFLE, mas um valor pequeno a= 0,05 foi simulado. 
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Figura 4.10. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo EEB, 
com F;=1,1 ecaso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura 4.11. Índices de sensibilidade das variáveis aleatórias à, para material CONC, modelo MDP, 


com F;=1,1 e caso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura 4.12. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MFLE, 
como F,=1,1 caso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 


As variáveis de resistência/oferta (G, e E) são identificadas por valores negativos 
de à; enquanto as variáveis de solicitação/demanda (P, a e C) por valores positivos dos 
cossenos diretores. Os gráficos indicam que a energia de fratura G, e o carregamento P, 
são as variáveis que mais influenciam na determinação do índice B para todos os valores 
de a e D em todos os modelos. É também evidente que todas as variáveis têm a mesma 
importância no cálculo de 8 para todos os tamanhos de viga (D). Uma exceção a esta 
tendência, e que o tamanho da trinca a, é ligeiramente mais importante para vigas mais 
espessas. Como esperado, o tamanho influencia significantemente a medida que a cresce, 
já que a= a,/D. É importante esclarecer que a, tem influência mesmo para a, =/D=0, já que 
o valor nulo para a indica somente o valor médio ua=0, mas com uma assumida variação 
(CV=0.10) como indica a Tabela 3.4. 

Finalmente, os valores de à; para C quantificam a maior influência do tamanho 
equivalente da ZP para o MDP. Isto se deve principalmente ao maior valor do CV quando 
aquela propriedade é calculada pelo MDP na Tabela 3.6. Nota-se que a influência de C 
é desprezível para o EEB devido ao baixo CV de C, segundo a Tabela 3.5 para material 
CONC. Todos os gráficos mostram que todas as variáveis consideradas como aleatórias têm 
uma importância no cálculo da confiabilidade e que nenhuma delas deve ser desprezada 
na análise de risco considerando a mecânica da fratura de vigas. 
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Similarmente, da Figura 4.13 á Figura 4.15 apresentam-se os valores dos cossenos 
diretores à; para o Caso Il. As mesmas tendências observadas para o caso | e inclusive as 
magnitudes de à; têm valores similares quando comparadas ao caso |, com um F =1,1 e 


para os F,=1,5 e 2,0, ilustrados no Anexo B. 
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Figura 4.13. Índices de sensibilidade das variáveis aleatórias à, para material CONC, modelo EEB, 
com F, = 1,1 e caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura 4.14. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MDP, 
com F,=1,1 e caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura 4.15. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MFLE, 


com F,=1,1 e caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 


A Tabela 4.1, Tabela 4.2 e Tabela 4.3 apresentam o resumo para o material CONC 


e para os fatores F, = 1,1;1,5 e 2,0, as tabelas com um critério de um 15% se analisa para 


cada parâmetro que influenciam no cálculo do índice B, sé tem incremento”, diminuiçãol, 


permanece constante o não tem influência com o incremento de a e D. 


Tabela 4.1. Influência dos índices de sensibilidade à; com F, = 1,1 para CONC e caso Il 


Parâmetro EEB MDP MFLE 
G,(N/m) aD! aD? a JDL 
E (Nim?) - - - 
C (mm) - aJDL - 
P(N) aD! aD? a JDT 
a, (m) a TDT a 1D7T atD7T 
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Tabela 4.2. Influência dos índices de sensibilidade à; com F, = 1,5 para CONC e caso Il 


Parâmetro EEB MDP MFLE 
G, (Nim) a DL a D*7 aJD? 
E (Nim?) - - - 
C (mm) - aJDL - 
PAN) aD! atD alD 
a, (m) a TDT a 1D7T a TDT 


Tabela 4.8. Influência dos índices de sensibilidade à; com F, = 2,0 para CONC e caso II. 


Parâmetro EEB MDP MFLE 
G, (Nim) aD! aD! a JDT 
E (Nim?) - - - 
C (mm) - aJDL - 
P(N) a* D! a D*L alD 
a, (m) a tD? a 1D? a TD) 


4.2 CONCRETO RECICLADO (FRAP) 


Os resultados de confiabilidade e de predição determinística de resistência para 
diferentes tamanhos de viga FTP e de trinca, variando os a para o material FRAP com Fç=1,1, 


avaliado como o algoritmo (HL-RF). 


4.2.1 Predição determinística do efeito de escala da resistência em vigas FTP 
para FRAP 

A predição da (0,) para o FRAP com relação à geometria da viga para os três 
modelos estudados é realizada da mesma maneira que para o CONC com as mesmas 
equações de desempenho. A mostra a dependência de com o tamanho das vigas, De o 
tamanho relativo das trincas (wc, =a,/D) para um de F,=1,1. Os gráficos foram gerados de 
acordo com a análise usual determinística, em que são utilizadas as médias y para o FRAP, 
com os parâmetros da Tabela 3.2 à Tabela 3.6. 

Para EEB e MDP, apresenta-se o mesmo efeito de escala de diminuição da o, com 
o aumento do tamanho da viga, como ilustrado nas tendências das curvas da Figura 4.16 
nos itens (a) e (b). Os dois modelos utilizam a mesma equação de o,, portanto as curvas 
de predição são similares e as diferenças são atribuídas aos parâmetros G,e C, que são 
maiores quando calculados pela metodologia EEB. 

Em vigas de FRAP para pequenos valores de a=0,1 com (D=0,063m), os modelos 
EEB e MDP também apresentam falhas de predição devido ao tamanho de ZP. As predições 
da MFLE, na Figura 4.16 para o item c), são bem conservadoras quando comparadas 
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às preditas pelo EEB e MDP. Além disso, as curvas de predição de resistência o, para o 
material (FRAP) tem um comportamento similar as curvas de resistência apresentadas por 
o material (CONC) na Figura 4.1, com a única diferencia de que o material (FRAP) ilustro 
valores de resistência menores em comparação com o material (CONC). 
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Figura 4.16. Predição deterministica da resistência o, e efeito de escala para a viga de FTP (FRAP): (a) 


EEB; (b) MDP e (c) MFLE. 


As curvas da resistência nominal com o aumento de tamanho em escala logarítmica 
são ilustradas na Figura 4.17. Assim como observado para o CONC, as curvas comportam- 


se como o efeito de escala clássico conforme Figura 2.5. 
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Figura 4.17. Predição deterministica da resistência o, e efeito de escala para a viga de FTP em escala 
logaritmica (FRAP): (a) EEB e (b) MDP. 


Com relação a comparação das resistências entre as vigas de CONC (Figura 4.1 
e Figura 4.2) e FRAP (Figura 4.16 e Figura 4.17). Como esperado, o efeito de escala é 
descrito pelo EEB e MDP através da diminuição da o, com o aumento do tamanho da 
espessura viga, como ilustram as tendências das curvas. As análises para estes materiais 
apresentam curvas com o mesmo comportamento e com as mesmas limitações para valores 
pequenos de D nos quais tem-se dificuldade com a convergência quando D aproxima-se 
à ZP, e as únicas diferenças estão nos valores dos parâmetros utilizados que são distintos 
para cada material. É importante notar que apesar dos materiais terem a mesma tendência, 
encontrarem-se valores de resistência maiores no CONC que para o FRAP, devido as 


características dos parâmetros dos materiais. 


4.2.2 Predição do índice de confiabilidade p e efeito de escala da resistência 
em vigas FTP 

Os índices B, para um foram determinados com o mesmo algoritmo e seguindo as 
mesmas equações de desempenho que no CONC, a Figura 4.18 apresenta os resultados 
de B para diferentes geometrias de viga e tamanhos relativo de trincas («,= a/D)e F,=1,1, 
como indicado na figura, as análises foram realizadas para as propriedades extraídas dos 
diferentes modelos EEB, MDP e MFLE para o caso Il de FRAP, que definem as variáveis 
aleatórias com uma distribuição lognormal (ver Tabela 3.2 á Tabela 3.6). O material FRAP 
analisou-se com a distribuição lognormal porque é uma de as mais usadas para fazer 
caracterização de parâmetros de laboratório, e porque está distribuição está limitada só 


para valores positivos. 
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Os índices de confiabilidade calculados para o EEB, MDP e MFLE demostram que 
devido ao efeito de escala, os índices B decrescem com o aumento da espessura da viga. 
Observa-se que o EEB tem um menor efeito de escala para B, como ilustram as curvas 
quase constantes. Os valores crescentes do tamanho da trinca relativos à altura, dados 
pela variável a, diminuem a confiabilidade e aumentam, desta forma, a probabilidade de 
falha. 
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Figura 4.18. Predição do índice de confiabilidade B para material FRAP, com Fç = 1,1 e Caso Il: (a) 
EEB; (b) MDP e (c) MFLE. 


É relevante notar que MFLE apresenta maiores limitações na convergência para 
a>0 e a=0,1 para valores pequenos de D, onde a espessura aproxima-se á ZP, como pode 
evidenciar-se na Figura 4.18 para o item c. 

Os índices de B, tanto para o CONC quanto para FRAP, seguiram análises com o 
mesmo algoritmo e com as mesmas equações de desempenho, variando os parâmetros 
característicos de cada material e apresentaram curvas com trajetórias similares e com 
incidência do efeito de escala, onde os valores de B decrescem com o aumento da 
espessura da viga para os dois materiais. O modelo de EEB para os dois materiais ilustra 
curvas constantes que são conservadoras em comparação com os outros dois modelos 
MDP e MFLE. Para todos os modelos os maiores índices B forem identificados para as 
menores relações de a, e para os menores valores da espessura da viga D. 

Para uma comparação com os resultados para o CONC, a Figura 4.19 apresenta os 
valores de B para F,=1,1 e a =0,25. Note que as curvas do FRAP têm um comportamento 
similar as curvas do CONC, mas com índices de confiabilidade menores para todos os 
modelos (EEB, MDP e MFLE). 
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Figura 4.19. Índices de confiabilidade CONC e FRAP caso Il com a=0,25 e para Fç=1,1;1,5€2,0: (a) 
EEB; (b) MDP e (c) MFLE 


A Figura 4.20 em 3D para o material FRAP com F, = 1,1 e caso Il apresenta os 
índices B em função de D e a, estabelecendo comparação entre EEB e MDP para os quais 
se ilustra uma incidência do efeito de escala de Bazant de uma diminuição dos índices B, 
com o incremento do de De a. 

Com índices B maiores para o MDP em comparação com EEB, um comportamento 
similar ao do material CONC para o mesmo valor de F, = 1,1. 
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Figura 4.20. Índices de confiabilidade B com Fs= 1,1 para os modelos EEB e MDP, material FRAP 
(caso II) 


4.2.3 Índice de sensibilidade à, das variáveis na determinação da confiabilidade 


Os índices de sensibilidade para o concreto reciclado estão definidos pela Equação 
(4.4) e são apresentados na Figura 4.21 á Figura 4.23 para o Caso Il do FRAP. Os à; ou 
cossenos diretores apresentam uma tendência de valores negativos para as variáveis G,e 
E, e uma tendência de valores positivos para as variáveis P., a, e C. A energia de fratura G, 
e o carregamento P, são as variáveis que mais influenciam as análises de confiabilidade. 
É importante esclarecer que a tem influência mesmo para «= a0/D=0, já que o valor nulo 
para a indica somente o valor médio ua=0, mas com uma assumida variação (CV=0,10), 
como indica a Tabela 3.4. 

Finalmente, os valores de à; para C quantificam a maior influência da zona de 
processo para o MDP. Isto se deve, principalmente, ao maior valor do CV quando aquela 
propriedade é calculada pelo MDP (Tabela 3.6) para o concreto reciclado. Observa-se que 
a influência de C é desprezível para o EEB devido ao baixo COV de C, de acordo com a 
Tabela 3.5. Ainda alguns dos parâmetros analisados tenham maior influência o maiores 
índices de sensibilidade, todos são de grande importância e ninguém de eles deverá de 
ser desprezado. 
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Figura 4.21. Índices de sensibilidade das variáveis aleatórias à; para material FRAP, modelo EEB, com 
Fç=1,1 e caso Il: (a) 050; (b) a=0.1; (c) a=0.25 e (d) a=0.4. 
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Figura 4.22. Índices de sensibilidade das variáveis aleatórias à; para material FRAP, modelo MDP, com 
Fo =1,1 caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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mD1=0,063 mD2=0,15 = D3=0,25 = D4=1,0 E D1=0,063 E D2=0,15 = D3=0,25 E D4=1,0 


0,80 0,80 

0,60 0,60 

0,40 0,40 + 
0,20 2 0,20 4 
1 0,00 E] a pri 0,00 = — 
8,020 “8 -0,20 
0,40 5 0,40 

-0,60 0,60 

-0,80 0,80 

-1,00 - 1,00 

Gf(N/m) E(N/m2) Ps(N) a(m) Gf(N/m) E(N/m2) Ps(N) alm) 
(a) a=0 (b) a=0,1 

0,80 E D1=0,063 m D2=0,15 = D3=0,25 E D4=1,0 0,80 + E D1=0,063 E D2=0,15 E D3=0,25 E D4-1,0 

0,60 0,60 

0,40 0,40 + 
s 0,00 = E E Em = 4 0,00 E 
ç -0,20 E 0,20 
w -0,40 n -0,40 

0,60 -0,60 

0,80 0,80 

1,00 1,00 

GfIN/m) E(N/m2) Ps(N) a(m) GfN/m) E(N/m2) Ps(N) a(m) 
(c) a=0,25 (d) a=0,4 


Figura 4.23. Índices de sensibilidade das variáveis aleatórias à; para material FRAP, modelo MFLE, 
com F,=1,1 e caso Il: (a) a=0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 


A Tabela 4.4 apresentam o resumo para o material FRAP e para os fatores F = 1,1, 
com um critério de um 15% se analisa para cada parâmetro que influenciam no cálculo do 
índice B, sé tem incremento”, diminuição!, permanece constante o não tem influência, 
com o incremento de a e D. 


Tabela 4.4. Influência dos índices de sensibilidade à; com F,= 1,1 para FRAP e caso II. 


Parâmetro EEB MDP MFLE 
G, (Nim) aD! at DT alD 
E (Nim?) - - - 

C (mm) - aJDL - 
PAN) - a D*1 alD 
a, (m) atD? a 1D? a TDT 
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ANÁLISE PROBABILÍSTICA DA RESISTÊNCIA DE VIGAS PELO 
MÉTODO DE MONTE CARLO (MMC) 


Neste capítulo as funções de densidade de probabilidade (FDP) da resistência de 
vigas em FTP são obtidas usando a técnica de simulação do método de Monte Carlo (MMC). 
As análises com o MMC utilizam o metodologia descrita na Seção 2.7 com os parâmetros 
de laboratório ensaiados para os concretos CONC e FRAP. As simulações foram realizadas 
seguindo a mesma equação de desempenho para o EEB e para MDP dada pela Equação 
(2.25) de 0,, para um número de simulações N = 50000. 

No MDP, os tipos de distribuições que seguem os parâmetros G,e CTOD,. foram 
determinadas por meio de simulações de números aleatórios e um teste de aderência a 
partir dos dados de laboratório do ensaio do MDP. O parâmetro Cfoi calculado indiretamente 
dos resultados daqueles parâmetros por meio das equações descritas na Seção 2.3.2, 
especificamente a Equação (2.43). No modelo EEB, as distribuições de C e G, foram 
determinadas por meio da simulação de números aleatórios, para as equações da Seção 
2.3.1, e posteriormente com um teste de aderência. 

O método de Monte Carlo para os EEB e para MDP, foi realizado por meio da 
variação geométrica da viga com a (0,05; 0,1; 0,40) e D (0,063; 0,15; 0,50; 1,0 m), sendo 
«q = a/D. 

Para parâmetros de fratura determinados em laboratório e para a simulação de 
resistência da viga FTP, forem feitos testes de aderência com o programa de matlab 
gerando números aleatórios com a função random para N=50000 simulações para obter a 
FDP, e depois por meio da ferramenta dfittool foi possível conhecer o tipo de distribuição 
que melhor ajustava-se a FDP de cada parâmetro analisado. É importante mencionar que 
os testes de aderência estão no anexo D. 


5.1 PARÂMETROS DE FRATURA DETERMINADOS EM LABORATÓRIO 


A Tabla 5.1 apresenta os dados dos parâmetros de laboratório que forem utilizadas 
para gerar a simulação do MMC para determinar as características probabilísticas e tipos de 
distribuições dos parâmetros Ce G,, para o EEB, e C, G,e CTOD,, para o MDP. Note que a 
4, O e tipo de distribuição dos parâmetros da Tabela 5.1 foram determinados baseados nos 
resultados dos ensaios das vigas como descrito na Seção 3.1. e na Seção 3.2. para CONC. 
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Tabela 5.1. Parâmetros estatísticos obtidos de laboratório para as simulações do MMC 


Variáveis V fo) Distribuição 
g(a,) 48,18 3,60 Normal 
gta) 207,04 14,02 Normal 
BN) 3,14E-12 0,49E-12 t-student 
A (Nim) 1,97E 0,29E-11 t-student 
E (N'm?) 34,82E+09 0,07E+09 lognormal 


A partir dos dados estatísticos da Tabela 5.1, as funções densidade de probabilidade 
(FDP) para a energia de fratura G, para os modelos EEB e MDP são apresentadas na 
Figura 5.1. As densidades foram determinadas por meio de N = 50000 simulações do MMC 
para os modelos EEB e MDP segundo a Equação (2.21) para o método de EEB e a partir de 
um teste de aderência dos dados obtidos do laboratório para G,no MDP respectivamente. 
Os parâmetros estatísticos das FDP são apresentados na Tabela 5.2 para os dois modelos. 
O tipo de distribuição indicado foi determinado pela distribuição que melhor se adaptou a 
curva da distribuição simulada. 


045 T 


EEB 
===MDP 


G, [Nm] 


Figura 5.1. Função densidade de probabilidade (FDP) do parâmetro G, para modelos: EEB (distribuição 
t-student) e MDP (distribuição normal). 


Embora os dois modelos tenham PDF de tipos diferentes: distribuição t-student 
para G, segundo o EEB, e distribuição normal quando calculado segundo as equações do 
MDP, ambas as curvas são similares em sua forma. A distribuição t-student faz parte das 
distribuições de probabilidade contínuas, onde as curvas de densidade tem uma distribuição 
com forma de sino com uma tendência similar à distribuição normal padrão. 
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Tabela 5.2. Parâmetros das FDP de G, (N/m). 


Parâmetros EEB MDP 
DZ) 70,392 61,740 
o 10,357 9,740 
Cv 0,147 0,150 
v, 3,457 - 

Distribuição t-student normal 


Similarmente, as funções densidade de probabilidade (FDP) para o tamanho 
equivalente da zona de processo, C, para os modelos EEB e MDP são apresentadas 
na Figura 5.2. Estas densidades também foram determinadas por meio de N = 50000 
simulações do MMC para os modelos EEB e MDP segundo as Equações (2.22) e (2.43), 
respectivamente. Os parâmetros estatísticos das FDP são apresentados na Tabela 5.3 
para os dois modelos. O tipo de distribuição indicado foi determinado pela distribuição que 
melhor se adaptou a curva da distribuição simulada. 


50 T T T T T 


—— EEB 
===MDP 


j l 
0 0,02 0,04 0.05 0,08 0,1 0.12 
€ [m] 


Figura 5.2. Função densidade de probabilidade (FDP) do parâmetro C para os modelos: EEB 
(distribuição t-student) e MDP (distribuição Weibull). 


Os dois modelos induzem PDF de tipo e forma distintas: distribuição t-student para 
C calculado pelo EEB, e distribuição Weibull quando calculado segundo as equações do 
MDP. É importante observar também a diferença marcante entre as médias para as duas 
distribuições, com um valor significantemente maior para a média de C pelo EEB. A maior 
dispersão da FDP gerada pelo MDP também pode ser observada na curva e pelos valores 
maiores dos parâmetros ce CV. 
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Tabela 5.3. Parâmetros das FDP de € (mm). 


parâmetros EEB MDP 
DZ) 37,000 21,000 
fo] 8,000 14,000 
Cv 0,216 0,667 
A 4,951 - 

A - 23,000 
Y = 1,397 
Distribuição t-student Weibull 


A principal diferença entre as metodologia para a determinação de C determinado 


pelo MDP, é que este é calculado através da propriedade de fratura CTOD, de acordo 


com a Equação (2.39). Desta forma, a Figura 5.3 apresenta a FDP com os valores dos 


parâmetros da curva descritos na Tabela 5.4. 


50 


45H 


40 


35 - 


30 — 


FDP 
x 
T 


20 — 


SER 


1 j 
0,01 0.02 


l j j 1 
0,03 0,04 0,05 0,06 0,07 0,08 0.09 


CTOD [m] 


Figura 5.3. Função densidade de probabilidade (FDP) do parâmetro CTOD,, para o MDP (distribuição 


normal). 


Tabela 5.4. Parâmetros da FDP de CTOD,, (m) do MDP. 


Parâmetros MDP 
DZ; 0,016 
o 0,006 
Cv 0,375 
Distribuição Normal 


Análise de confiabilidade em vigas FTP 


71 


5.2 SIMULAÇÃO DA RESISTÊNCIA (0) DA VIGA FTP 


Com base nos parâmetros estatísticos para as propriedades de fratura de cada 
modelo (EEB e MDP) da seção anterior, simulações da resistência, o,, da viga FTP foi 
realizada por meio do MMC com 50000 combinações aleatórias simuladas para a predição 
da resistência da viga (ver Figura 2.4) segundo a Equação (2.25). Diferentes combinações 
de geometria foram consideradas para a viga FTP com a (0,05; 0,1; 0,40) e D (0,063; 
0,15; 0,50; 1,0 m). Os parâmetros estatísticos para a geometria da viga adotados para as 
simulações são mostrados na Tabela 5.5. 


Tabela 5.5. Parâmetros geometricos para a viga FTP a serem usados para simulações do MMC 


Variáveis V fo) Distribuição 
a (m) (0,05; 0,1; 0,40) 0,10 Lognormal 
Dm) [0,063;0,15;0,5;1,0] -— Determinístico 
t(m) 0,08 -— Determinístico 


As curvas das FDP adaptadas aos dados simulados pelo MMC são apresentadas 
da Figura 5.4 à Figura 5.7 para tamanhos de viga D de 0,063 m, 0,15 m, 0,50 m e 1,0 
m, respectivamente. Para cada figura de um determinado tamanho, as FDP para os 
valores médios de a, como indicados na Tabela 5.5, são plotadas para ambos os modelos 
analisados: EEB e MDP. Assim pode-se comparar diretamente as curvas geradas para 
diferentes tamanhos relativos de trinca (a) e modelos de fratura em um determinado 
tamanho D. 

Os parâmetros estatísticos das FDP ilustradas nas figuras são descritos 
quantitativamente da Tabela 5.6 à Tabela 5.13 para todos os tamanhos analisados e os 
dois modelos. As tabelas também indicam o tipo de distribuição que melhor se adaptou aos 
dados simulados pelo MMC. 
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Figura 5.4. Função densidade de probabilidade (FDP) de o, para EEB e MDP com D = 0,063 m. 


Tabela 5.6. Parâmetros da FDP de o, (N/m?) para EEB com D=0,063 m. 


a=0,05 a=0,1 a=0,4 
DZ) 2,89E +06 2,89E+06 1,49E+06 
[oi 0,38E+06 0,35E+06 0,22E+06 
CV 0,131 0,121 0,148 
V, 2,896 2,923 4,154 
Distribuição t- Student t- Student t- Student 


Tabela 5.7. Parâmetros da FDP de o, (N/m?) para MDP com D = 0,063 m. 


a=0,05 a=0,1 a=0,4 

3,35E+06 3,79E+06 2,01E+06 

1,02E+06 1,10E+06 0,49E+06 
Cv 0,304 0,290 0,243 
uy - 15,108 14,487 
0, - 0,283 0,239 
k, 0,090 - ad 
Distribuição VEG. Lognormal Lognormal 
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= 
—— EEB 
— —»MDP 


FDP 


o, [Nim] x10º 


Figura 5.5. Função densidade de probabilidade (FDP) de o, para EEB e MDP com D= 0,15 m. 


Tabela 5.8. Parâmetros da FDP de o, (N/m2) para EEB com D=0,15 m. 


a=0,05 a=0,1 a=0,4 
7) 2,72E+06 2,60E+06 1,26E+06 
fo] 0,32E+06 0,29E+06 0,18E+06 
Cv 0,117 0,111 0,143 
v, 3.005 3,380 5,104 
Distribuição t- student t- student t- student 


Tabela 5.9. Parâmetros da FDP de o, (N/m?) para MDP com D= 0,15 m. 


a=0,05 a=0,1 a=0,4 
V 3,51E+06 3,11E+06 1,55E+06 
fo] 0,96E+06 0,66E+06 0,29E+06 
Cv 0,27 0,212 0,187 
uy, 15,035 14,929 14,236 
0, 0,269 0,209 0,184 
Distribuição Lognormal. Lognormal Lognormal 


Análise de confiabilidade em vigas FTP 80 


2 25 
o, [Nim] 


Figura 5.6. Função densidade de probabilidade (FDP) de o, para EEB e MDP com D =0,5 m 


Tabela 5.10. Parâmetros da FDP de o, (N/m?) para EEB com D= 0,5 m. 


a=0,05 a=0,1 a=0,4 
DZ) 2,26E+06 1,96E+06 0,85E+06 
fo] 0,23E+06 0,18E+06 0,11E+06 
Cv 0,102 0,092 0,129 
v, 3,600 3,864 5,541 
Distribuição t- student t- student t- student 
Tabela 5.11. Parâmetros da FDP de o, (N/m?) para MDP com D=0,5 m. 
a=0,05 a=0,1 a=0,4 
DZ) 2,56E+06 2,10E+06 0,96E+06 
fo] 0,43E+06 0,26E+06 0,13E+06 
CV 0,168 0,124 0,135 
Distribuição Normal Normal Normal 
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t5 


o, [Nim] 


Figura 5.7. Função densidade de probabilidade (FDP) de o para EEB e MDP com D=1,0 


Tabela 5.12. Parâmetros da FDP de o, (N/m?) para EEB com D=1,0 m. 


a =0,05 a=0,1 a=0,4 
24 1,87E+06 1,53E+06 0,64E+06 
fo] 0,17E+06 0,13E+06 0,08E+06 
[ONA 0,090 0,085 0,125 
v, 3,835 3,799 5,448 
Distribuição t- Student t- Student t- Student 


Tabela 5.13. Parâmetros da FDP de o, (N/m?) para MDP com D = 1,0 m. 


a =0,05 a=0,1 a=0,4 
DZ) 2,01E+06 1,58E+06 0,70E+06 
fo] 0,25E+06 0,16E+06 0,09E+06 
Cv 0,124 0,101 0,129 
Distribuição Normal 


Normal 


Normal 


Análise de confiabilidade em vigas FTP 


82 


Ao comparar as curvas geradas para as diferentes geometrias nos modelos de EEB 
e MDP, pode-se observar que a resistência o, para o EEB seguem predominantemente 
distribuições t-student para todos os tamanhos relativos de trinca e tamanhos de viga. Por 
outro lado, o MDP apresentou distribuições variadas, dentre elas valor extremo generalizado 
(VEG), lognormais, e normais. 

Analisando especificamente as PDF para o modelo EEB, a resistência o, segue 
distribuição t-student com o parâmetro grau de liberdade, v,, crescente com o aumento de a. 
É importante observar que a distribuição t-student é bem similar a distribuição normal, mas 
apresenta as caudas (tail da distribuição) mais alongados. Ademais, as PDF apresentam 
uma crescente dispersão com o aumento de D. 

Para as distribuições do modelo de MDP, a resistência o apresenta, na maioria dos 
casos (diferentes D e a) distribuições normais e lognormais. Por outro lado, a distribuição 
de VEG foi melhor adaptada para o casos em que o tamanho viga era menores (D=0,063 
m). Distribuições de valor extremo são comumente associadas a modelagem de problemas 
com probabilidades muito pequenas ou eventos raros. 
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CONCLUSÕES E SUGESTÕES PARA TRABALHOS FUTUROS 


6.1 CONCLUSÕES 


Esta pesquisa implementou com sucesso um algoritmo iterativo para o cálculo 
do índice B, para funções de estados limites não lineares. O algoritmo mostrou bom 
desempenho e convergência para todos os cenários analisados. Os resultados mostraram 
que a análise determinística do concreto utilizando parâmetros de fratura extraídos de 
ensaios com as metodologias EEB e MDP são similares para a predição do efeito de escala 
da resistência nominal de vigas elaboras em os materiais CONC e FRAP. Entretanto, as 
predições dos índices de confiabilidade considerando as incertezas dos parâmetros do 
material, geometria e carregamento foram bastante distintas para cada um dos modelos. 

As confiabilidades determinadas com parâmetros do EEB foram bem mais 
conservadoras (menor valor de B) para os diversos tamanhos de trinca a e espessuras de 
viga D quando comparados com as predições de confiabilidade a partir dos parâmetros de 
MDP. Além disso, este último mostrou uma grande variação de B para a variação de a e D 
considerados. Isto é explicado principalmente pelo mais alto CV do parâmetro C quando 
determinados pelo MDP. Este alto CV é atribuído a grande incerteza de C que é calculado 
indiretamente por várias equações não lineares e em medidas de inclinações de curvas de 
histeresis em ensaios do material. Enquanto que C, no modelo EEB, é calculado baseado 
em medidas de carregamentos máximos em vigas de tamanhos diferentes que tendem a 
ter menor variação em suas medidas. 

As análises apresentadas também mostram a existência de um efeito de escala 
na tendência do índice de confiabilidade, em que estruturas maiores (maior D) tem maior 
probabilidade de falha. Além do mais, a variabilidade do parâmetro C é a que mais afeta 
esse efeito de escala de B. Este trabalho também demonstrou quantitativamente que a 
escolha do modelo, e as metodologias experimentais associadas a ele, contribui de maneira 
significativa para a determinação da confiabilidade estrutural. Além disso, os resultados 
aqui mostram que a consideração da variabilidade estatísticas destes parâmetros do 
material associados as incertezas do carregamento e geometria da peça são essenciais 
para a adequada predição de falha estrutural. 

O conhecimento de os índices de confiabilidade, probabilidades de falha e incidência 
dos parâmetros do concreto para variação geometria da viga de FTP com fraturamento 
tipo | que é um tipo de falha muito comum, servirá aos calculistas como guia nos projetos 
de dimensionamento estrutural considerando a mecânica da fratura. Com a utilização de 
diferentes coeficientes de segurança (F.) na equação de estado limite para o EEB, MFLE 
e MDP foi possível concluir que para maiores fatores de segurança obtém-se maiores 
índices de confiabilidade é dizer que estão relacionado de maneira direita. Ao estabelecer 
comparações entre os dois materiais estudados, o CONC e o FRAP com a avalição 
do método de confiabilidade, determinou-se que os dos apresentarem gráficos com as 
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mesmas tendências e com a mesma incidência do efeito de escala, as únicas diferencias 


que encontraram-se estavam em relação que para o CONC os índices p e a resistência 


resultarem maiores que para o FRAP. 


As simulações do MMC por meio da geração de números aleatórios permitiu 


conhecer as distribuições que apresentavam as variáveis que usadas nos modelos EEB e 


MDP descrevendo o comportamento probabilístico da resistência para o material CONC. 


As FDP de do EEB apresentarem uma predominância da distribuição t-Student, e no caso 


do MDP, predominou a distribuição normal e lognormal. 


6.2 SUGESTÕES PARA TRABALHOS FUTUROS 


Como sugestões para trabalhos futuros, pode-se citar: 


determinação do índice de confiabilidade para outras estruturas e condições de 
carregamentos: placas planas, punção em lajes e outras condições; 


estudo o efeito do CV de G,e C no valor do índice de confiabilidade para vários 
tamanhos; 


formulação de outras equações de estado limite para fratura de estruturas de 
concreto; 


consideração outros tipos de distribuição para as variáveis geométricas da viga 
além das distribuições normais e lognormais. 
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ANEXOS 


A. ÍNDICES DE CONFIABILIDADE. 


Analise de confiabilidade para F, = 1,5 e F, = 2,0 concreto virgem, índices B para 
distribuições normais e lognormais caso | e caso Il, para material CONC. 
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Figura A.1. Predição do índice de confiabilidade B para material CONC, com F, = 1,5 e caso |: (a) EEB; 
(b) MDP e (c) MFLE. 
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Figura A.2. Predição do índice de confiabilidade B para material CONC, com F, = 1,5 e caso Il: (a) EEB; 
(b) MDP e (c) MFLE 
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Figura A.4. Predição do índice de confiabilidade B para material CONC, com F, = 2,0 e caso Il: (a) EEB; 
(b) MDP e (c) MFLE 
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B. ÍNDICES DE SENSIBILIDADE 


Analise de confiabilidade para F, = 1,5 e F, = 2,0, cossenos &; para distribuições 


normais e lognormais caso | e caso Il, para material CONC. 
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Figura B.1. Índices de sensibilidade das variáveis aleatórias à, para material CONC, modelo EEB, com 


Fo=1,5e caso I: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.2. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MDP, com 
Fç=1,50caso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.3. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MFLE, 
com F.=1,5e caso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.4. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo EEB, com 
F;=1,5 e caso Il: (a) a50; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.5. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MDP, com 
Fç=1,5 e caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.6. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MFLE, 
com F, = 1,50 caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.7. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo EEB, com 
F.=2,0ecaso: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 


Anexos 101 


E D1=0,063 m D2=0,15 = D3-0,25 m D4=1,0 


0,80 4 080 mD1=0,063 8 D2=0,15 = D3=0,25 mD4=1,0 
0,60 + 0,60 
0,40 + 9,40 
a 
a 0,20 + “2 0,20 
- Ss 
S om o .? Ro o. 
& 3 
F -0,20 + E020 
o) E) 
w -D,40 4 tn 0,40 
-0,60 + -0,60 
-0,80 + -0,80 
-1,00 + -1,00 
Gf(N/m) E(N/m2) Clm) Ps(N) a(m) Gf(N/m) E(N/m2) clm) Ps(N) am) 
(a) 450,0 (b) a=0,1 
ED1-0,063 8 D2=0,15 mD3-0,25 8D4-1,0 oo = mD1-0,063 1 D2=0,15 =D3-0,25 1 D4-1,0 
0,80 
0,60 4 0,60 
0,40 0,40 + 


Eq Ea di sl 
je 
É 900 me O a 0,00 


Ea 
zo —— au ht 
3 «5 
No) s=2-0,20 + 
s$-020 É 
&-0,40 -0,40 + 
-0,60 -0,60 + 


GF(N/m) E(N/m2) clm) Ps(N) alm) Gf(N/m) E(N/m2) clm) Ps(N) a(m) 


(0) d=0,25 (d) a=0,4 


Figura B.8. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MDP, com 
Fo=2,0 e caso I: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.9. Índices de sensibilidade das variáveis aleatórias à; para material CONC, modelo MFLE, 
com F.=2,0 e caso |: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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Figura B.10. Índices de sensibilidade das variáveis aleatórias &; para material CONC, modelo EEB, 
com F,=2,0 e caso Il: (a) a>0; (b) a=0.1; (c) a=0.25 e (d) a=0.4. 
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Figura B.11. Índices de sensibilidade das variáveis aleatórias à, para material CONC, modelo MDP, 
com F.=2,0 e caso Il: (a) a>0; (b) a=0,1; (c) a=0,25 e (d) a=0.4. 
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Figura B.12. Índices de sensibilidade das variáveis aleatórias É; para material CONC, modelo MFLE, 
com FS = 2,0 e caso Il: (a) a=0; (b) a=0,1; (c) a=0,25 e (d) a=0,4. 
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C. FATORES DE SEGURANÇA 
Tabelas para os F = 1,1; 1,5 e 2,0 para CONC e F = 1,1 para FRAP 


Tabela C.1. Índices de confiabilidade B com Fo =1,1 para CONC e FRAP e caso Il. 


Fs=1,1 CONC FRAP 

EEB MDP MFLE EEB MDP MFLE 
D(m) B a B a B a B a B a B a 
0,063 0,6671 0,00 1,0816 0,00 N.c 0,00 0,3560 0,00 0,8444 0,00 N.c 0,00 
0,15 0,6673 1,0750 “ 0,3561 0,8386 - 
0,25 0,6676 1,0675 “ 0,3562 0,8320 - 
1,00 0,6688 1,0129 - 0,3570 0,7834 a 
0,063 0,6686 0,10 0,9881 0,10 7,1029 0,10 0,3566 0,10 0,7604 0,10 7,1029 0,10 
0,15 0,6695 0,8748 4,7697 0,3579 0,6515 4,7697 
0,25 0,6689 0,7736 3,6397 0,3587 0,5459 3,6397 
1,00 0,6617 0,4634 1,6705 0,3597 0,1697 1,6705 
0,063 0,6611 0,25 0,8913 0,25 4,924 0,25 0,36344 0,25 0,6848 0,25 4,924 0,25 
0,15 0,6479 0,7139 3,1545 0,3605 0,5056 3,1545 
0,25 0,6409 0,5951 2,3686 0,3588 0,3665 2,3686 
1,00 0,6298 0,3482 1,1582 0,3557 0,0339 1,1582 
0,063 0,6144 0,40 0,8045 0,40 3/7104 0,40 0,3613 0,40 0,6515 0,40 3,7104 0,40 
0,15 0,5788 0,6102 2,4469 0,3498 0,4556 2,4469 
0,25 0,5649 0,5006 1,8618 0,3441 0,3192 1,8618 
1,00 0,5502 0,2982 0,9491 0,3367 0,0241 0,9491 
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Tabela C.2. Índices de confiabilidade B com F, = 1,5 para CONC e caso II. 


F=1,5 CONC 

EEB MDP MFLE 
D(m) B a B a B a 
0,063 2,9299 0,00 2,9797 0,00 N.c 0,00 
0,15 2,9302 2,9813 . 
0,25 2,9304 2,9829 “ 
1,00 2,9304 2,9902 “ 
0,063 2,9262 0,10 2,9927 0,10 9,2525 0,10 
0,15 2,9228 2,9803 6,9221 
0,25 2,9141 2,9446 5,7932 
1,00 2,8654 2,7154 3,8254 
0,063 2,8383 0,25 2,9085 0,25 6,8136 0,25 
0,15 2,7697 2,7865 5,1199 
0,25 2,7326 2,6863 4,3570 
1,00 2,6745 2,4498 3,1748 
0,063 2,5498 0,40 2,6353 0,40 4,9584 0,4 
0,15 2,3675 2,3636 3,9004 
0,25 2,2966 2,2295 3,3907 
1,00 2,2240 2,0183 2,5751 


Tabela C.3. Índices de confiabilidade B com Fo = 2,0 para CONC e caso Il 


F=20 CONC 
EEB MDP MFLE 

D(m) B a B a B a 

0,063 5,0287 0,00 4,7391 0,00 N.c 0,00 

0,15 5,0291 4,7469 N.c 

0,25 5,0293 4,7556 - 

1,00 5,0281 4,8117 - 

0,063 5,0309 0,10 4,8360 0,10 N.C 0,10 

0,15 5,0193 4,9130 - 

0,25 5,0000 4,9430 - 

1,00 4,9073 4,8027 - 

0,063 4,8490 0,25 4,7719 0,25 8,4383 0,25 

0,15 4,7155 4,6928 6,8678 

0,25 4,6422 4,5974 6,1404 

1,00 4,5303 4,3540 4,9995 

0,063 4,2799 0,40 4,3303 0,40 5,8795 0,4 

0,15 3,8998 3,8844 5,0332 

0,25 3,7555 3,6859 4,6107 

1,00 3,6196 3,4500 3,9143 
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D. PARÂMETROS SIMULADOS POR MONTECARLO 


Testes de aderência realizados para os parâmetros C, CTMOD e G, e para as os 
diferentes combinações geométricas de a (0,05; 0,1; 0,40) e D (0,063; 0,15; 0,50; 1,0 m) 


realizados para as distribuições geradas pelo MMC para material CONC. 


FDP 


Clml 


(a) EEB 


Clm] 


(b) MDP 


Figura D.1. Função densidade de probabilidade (FDP) do parâmetro C para os modelos: a) EEB 
(distribuição t-student) e b) MDP (distribuição Weibull). 
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G, [N/m] 


(a) EEB 


50 


60 7 
G, [Nim] 


(b) MDP 


100 


Figura D.2. Função densidade de probabilidade (FDP) do parâmetro G, para modelos: a) EEB 
(distribuição t-student) e b) MDP (distribuição normal). 
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CTOD [ml 


Figura D.3. Função densidade de probabilidade (FDP) do parâmetro CTOD,, para o MDP (distribuição 
normal). 
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Figura D.4. Parâmetro com D=0,063m com a =0,05, a=0,1 e a=0,40: para EEB a),e) e f) distribuição 
T-Student e para MDP: c) valor extremo generalizado d) e f) Lognormal. 
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Figura D.5. Parâmetro o, para D=0,15 m com a =0,05, a=0,1 e a=0,40: para EEB a), e) e f) distribuição 
t-student e para MDP: c), d) e f) Lognormal. 
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Figura D.6. Parâmetro o, para D =0,5 m com a =0,05, a=0,1 e a=0,40: para EEB a), e) e f) distribuição 


t-student e para MDP: c), d) e f) Normal. 
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Figura D.7. Parâmetro o, para D =1,0 m com a =0,05, a=0,1 e a=0,40: para EEB a), e) e f) distribuição 
T-Student e para MDP: c), d) e f) Normal. 
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